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REXOR II ROWRCRAFT SIMULAnON MODEL* 

Volume II CcHsputer Implenentation 
J. S. Reaser and P. H. Kretsinger 
Lockheed-Califbrnia Company 
SUMMARY 

This report describes a generalized format rotorcraft nonlinear 
simulation called REXOR II. Ihe program (models single main rotor vehicles 
with up to seven main rotor blades. Wings, two horizontal tail planes, 
and auxiliary thrustors may be included to model a variety of conrpoand 
helicopter configurations. 

Program output is primarily in the form of machine plotted time his- 
tories specified ‘from a signal list. This list is, in turn, user selected 
from a set of computation variables used by the program. 


LIST OF SYMBOLS 


The symbols used in the REXOR II equations are quite numerous. In order to 
keep the catalog of symbols to manageable proportions a list of basic symbols 
is given, followed subscripts, superscripts, and postscripts. Nonconform- 
ing cases of usage together with c<mplicated or obscure subscripting are fully 
annotated in the basic list. 

SYMBOLS 


a. 


1 

[Aj 


1 , 2,3 


In 


arbitrary vector 
speed of sound, m/s 

2 

acceleration vector, m/s (ft/s ) 

longitudinal component of blade first harmonic flapping, rad 
generalized mass element matrix 
modal variables 

generalized displacement of n th blade, first mode 


*The contract reseeurch effort which has lead to the results in this report 
was financially supported by USARTL (AVRADCOM) Structures Laboratory, 



2n 

3n 

13 


b 


B 



c 


Cc] 


c 


D 


"M 


C 

P 

^X,Y,Z 


^ 1 , 2,3 

C(k) 

d 

dr 

dv; 

d/dt 

(d/e)o 

(d/e)^ 


e 


El 

f 

iMR 


generalized displacement of nth blade t second mode 
generalized displacment of nth blade* third mode 
cosine conponent of blade first harmonic cyclic, rad 
niunber of main rotor blades; arbitrary vector 
dissipation function 

sine component of blade first harmonic cyclic, rad 

blade segment chord, m (ft) 

damping matrix • 

aerodynamic drag coefficient 

aerodynamic lift coefficient 

aerodynamic pitching moment coefficient 

pover coefficient 

thrust coefficient 

linear damping, N/m/s (ib/ft/s) 

rotary damping, N-m/rad/s (ft- lb /rad/s) 

blade bending to feathering couplings 

lift deficiency function 

infinitesimal increment 

increment in rotor, radius, m (ft) 

increment in time 

derivative with respect to time 

swashplate to feather gear ratio, zero collective 

swashplate to feather gear ratio slope with collective 

pitch horn effective crauik arm, m (ft) 

p p 

blade bending stiffness distribution, N-m (ib-ft ) 
ground effect factor for main rotor 


2 



F 

F 

x.y.z 

F 

F 

ePH 

g 

®X,Y,Z 

G 

{-} 

G 

GJ 





J 

J 

k 



factor; force, N (lb) 

force coB^&ents along X,Y,Z directions, N (lb) 
genersdized force about 0, ♦ axis 
feathering node generalized force 
gravity, m/s^ (ft/s^) 

gravity components along X,Y,Z directions 
gear ratio 

generalized force vector 

gyro angular acceleration partial product 

blade torsional stiffness, N-m (lb - ft ) 

= Xj^^, kg-m^ (slug-ft^) 

= Yj^^, kg-m^ (slug-ft^) 

= kg-m^ (slug-ft^) 

= Em^ (Y^^ + Z^^), kg-m^ (slug-ft^) 

= Em. (X,^ + Z.^), kg-m^ (slug-ft^) 

= Em^ (X^^ + y^^), kg-m^ (slug- ft®) 

= Em^ Xj^ Yj, kg-m® (slug- ft®) 

= Em^^ X^ Z^t kg-m® (slug- ft®) 

= Em^ Yj^ Zj, kg-m® (slug- ft®) 
unit vector 
unit vector 
advance ratio 

number of blade radial stations; reduced frequency, 
rad/s; vinit vector 

spring matrix 

blade spring matrix element 



K. 


ft. 


X.Y.Z 


IB 


"OB 


n 

“f 

“h 


m. 

1 


m, 


SP 
M 

[m] 

«rk 

^“Sc.Y.Z 

M 

V 

N 

P 


^iMR 


spring constants along X,Y,Z direction, N/m (Ib/ft) 

spring rates about 0, 4» axis, N-m/rad (ft-lb/rad) 

location inboard feather bearing, m (ft) 

location outboard feather bearing, m (ft) 

radial location of intersection of precone and feather 
axis, o (ft) 

rolling maaent, N-b (ft- lb) 

mass of element, kg (slugs) 

summed fuselage coordinate mass, kg (slugs) 

summed hub axis mass, kg (slugs) 

mass of i th particle or blade segment, kg (slugs) 

svashplate summed mass, kg (slugs) 

pitching moment, N-m (ft-lb); = Em^, kg (slugs); mach number 
generalized mass matrix 
generalized mass matrix element 
* Em^ kg-tn (slug-ft) 

= Em^ Y^, kg-m (slug-ft) 

= Em^ Z^, kg-m (slug-ft) 

moments about X,Y,Z axis, N-m (ft-lb) 
blade torsional moment, N-m/m (ft-lb/ft) 
number of system particles 

angular velocity about X axis, rad/s; particle 

main rotor pitch mcxaent inflow, m/s (ft/s) 

generalized coordinate; angular velocity about Y axis, 
rad/s 
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‘^iMR 

Q 


QLO^DS 



r 


R 

R. 


Z4.Z 

S 

®HA 

t 

I 


C,P.S,T 


'^iMR 

*'iTR 

X 


sain rotor roll moaent inflow, m/s (ft/s) 
generalized forcing function 

aerodynamic pressure times reference wing area, kg (lb) 
tota^. nonoain rotor aerodynamic loads matrix 
tall rotor torque, N-m (ft-lb) 

general vector; radius of curvature, ft; angular velocity 
about Z axis, rad/sec; notation for (X,Y,Z) 

static blade shape 

vector displacement of particle p in X,Y,Z axis system 
vector displacement of x,y,z origin in X,Y,Z system 
gyro diuaper coupling ratios 

Laplace variable, path of notion of particle p 
blade spline length along neutral axis locii, m (ft) 
time 

kinetic energy, N-m (ft-lb) 
transformation of coordinates matrix 
velocity in X direction, m/s (ft/s) 

potential energy function, N-m (ft-lb); strain energy, N-m (ft-lb) 

air velocity on blade element, m/s (ft/sec) 

velocity in Y direction, m/s (ft/sec) 

trajectory velocity 

velocity in Z direction, m/s (ft/sec) 

main rotor collective inflow, m/s (ft/sec) 

tail rotor collective inflow, m/s (ft/sec) 

motion in X direction, ra (ft); blade span location 
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X 





J 

y 


1.2,3 


OSA 



"ITD 


1.2,3 


'03L 


^OF 


OSP 


6 


‘FA 


e 


PHr. 


r 


0 


coordinate direction; axis; teflecti<m, a (ft); location, 
a (ft); crt^s product 

blade radial station of sweep and Jog, a (ft) 
trajectory path, a (ft) 
tail rotor l<aigitudinal force, a (lb) 
motion in T directi<m, a (ft) 

coordinate direction; axis; deflecti<m, a (ft); location, a (ft) 
tension torsion padi ouUx>ard «d aodal coefficients 

difference between Y directi^ locations of eg and neutral 
axis points of blade elcaait, a (ft) 

action in Z direction 

coordinate directi<ai; axis; d«>' fleet i(», a (ft); location, a (ft) 

relative svashplate veirtical displaceaent with respect to 
the hub, a (ft) 

tension>torsion pack outboard end aodal coefficients 

teetering rotor undersling, a (ft) 

hub set distance above fuselage set, a (ft) 

hub set distance above swasfaplate set, a (ft) 

angle of attack, rad 

angle of attack with hub set, rad 

sideslip angle, rad 

blade feathering angle, rad 

feathering/pitch-hom bending or dynaiaic torsion 
generalized coordinate displaceaent 

blade droop relative to precone angle, rad 
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'T 

6 

*3TR 

9c/ 

d 


c 

T 



U 


blade sweep angle, rad; dyiunic stall delay, s 

trajectory path angle with E set, rad 

liBit deflection, rad; freeplay, rad; snail increment 

tail rotor pitch - flap coupling 

downvash factor of ving on horizontal tail 

vector notati<»i of 4, d, 4 

rotation about T axis, rad 

collective blade angle, rad 

sileslip at blate el^aent, rad 

air dmsity, kg/n^, (slugs/ft^) 

time coastant, s; natural period, s 
feathering axis precone, rad 
rotation about X axis, rad 
fathering angle, rad 

feathering angle of blade element of n th blade, rad 

blade root reference feather angle, rad 

blade torsion, rad 

SUB of blade tvist and torsion, rad 

vake angle of aain rotor, rad, (deg) 

rotation about Z axis, rad; sideslip angle with hub set, rad 

control input axis rotation from svashplate, rad 

pitch lead angle, rad, (deg) 

trajectory path yaw with E set, rad 

main rotor apparent airflow an^le, rad 

rotational speed, rad/s; angular velocity, rad/s; 
natural frequency, rad/s 
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3 


partial derivative, derivation 


SUBSCRIPTS 

a 


/. 


b 

BEND 

BLE 

BIji 

C 


CG 

com 

DW 

DYH 


EIIG 


EST 


FA 

FB 


rn 




G 


GSIJ 


GFB 

G3P 

a 


arbitrary coordinate set a 

due to aerodynaaics 

arbitrary coordinate set b 

associated with .lade elastic bending 

blade eleiaent coordinate systea 

blade reference axis systea for the n^ blade 

associated with pilot control input, chordwise 

associated with center of gravity location 

corrective, correction 

referring to downwash 

referring to dynamic component 

earth aucis 

aissociated with powerplant - engine 
estimated 

fuselage axis; associated with blade feathering 
referring to blade feather axis 
associated with feedback 

associated with feathering of the n th blade 
due to friction 

referring to gyro or gyro coordinate system 
associated with gas generator section of powerplant 
associated with gyro control feedback 
gyro to swashplate connection 



GUB relating to ©rro ginbal unbalance 

H referring to hub or principal reference axis syster. 

HT associated with horizontal tail 

i referring to inflow, particle 

IB referring to inboard feather bearing location 

J spring matrix index 

jog associated with blade attachment joggle 

J associated with gyro end of feedback rod linkage 


Jn 

k 

LAG 

LD!IT 

m 

MR 

n 


KA 


HEW 

HO 

NR 

OB 


OLD 


P 

PH 

r 

R 


associated with feedback rod coming from the n th blade 

generalized mass index 

associated with lead-lag damper 

signifying limiting value 

blade CK>de index, spring matrix index 

associated with main rotor 

blauie number index 

referring to blade segment neutral axis 

newly determined value 

normal (to airflow) component 

pertaining to nonrotating value 

referring to outboard feather bearing location 

value from previous time step 

associated with propeller; perpend icula.* blade component 
referring to pitch horn 
generalized mass index 
referring to rotor axis system 


a 



REF 

R. '-l 
S 

SC 

s? 

s? 

c 

S. SP 
STEADY 
S’./ 

T 

TR 

TRIM 


US 


UHSTEADY 




WING 


X 


Y 

YA 

Z 

Z.A 


associated with blade feather reference value 

referring to control gyro feedback lever moment 

refe?'ring to blade spanwise velocity; general mode; static; 
structural; shaft 

referring to blade segment shear center 

referring to swashplate 

command to swashplate 

referring to swashplate limit stop 

steady component 

referring to blade sweep angle location 

associated with trajectory path relating to E axis; 
tangential blade component; blade torsion; blade twist 

associated with the tail rotor 

initial or trim value 

associated with blade twist {built in) 

relating to control gyro unbalance 

associated with unsteady component 

associated with vertical tail 

associated with the wing 

relating tc component in X direction 

relating to component in Y direction 

relating to aerodynamic component in Y direction 

relating to component in Z direction 

relating to aerodynamic component in Y direction 
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0 


1.2.3 

IS 

1/k c 
3/1* c 

®PHn 

e 

♦ 

SUPERSCRIPTS 

I 

T 

{-) 

(’) 

(•) 

(••) 

(- 1 ) 

(-) 

POSTSCRIPTS 

(i) 

(n) 


(nought) associated with collective value, coordinate axis 
value, with respect to principal reference axis, blade 
root summation 

with respect to blade modes 1, 2, or 3 
first harmonic component shaft axis feathering 
with respect to blade 1/U chord 
with respect to blade 3/^ chord 

associated with the feathering mode of the n th blade 
relating to compouent in the direction 
relating to component in the 0 direction 
relating to component in the ^ direction 


referring to inertial reference 

matrix transpose 

(bar) average quantity 

(prime) slope with respect to blade span 

(dot) time derivative of basic quantity 

(double dot) second time derivative 

matrix inverse 

vector quantity 

blade radial station index 
blade number index 
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1. INTRODUCTION 


The presentation of the REXOR II software in this volume follows the 
multilevel hierarchial design of the program modules. A skeletal chart 
showing the major computational paths is presented in Figure 1. 

The modules are discussed as follows: 

# Section 2 — MAIN — Level 1 Executive 

e Section 3 “■ FLY — Level 2 Executive 

e Section U — ACCEL — Level 3 Executive 

e Section 5 ~ GFORCE — Level U Executive 

e Section 6 — MASSSP — Level I 4 Executive 

The design of REXOR II at the executive level is beised on the classical 
technique input-process-output. The process can be characterized as a signal 
generator, with the output being a defined set of time-varying functions 
(signals) and the input 'lock being ^at is necessary to define the physical 
model and control the process. The characterization of REXOR II as a signed 
generator provides a mechanism for extensive output processing. The signal 
set generated can be identified as input to a variety of output processors. 
Exaiiq)les are: 

e Printed time histories 

e Graphic displays 

e Haurmonic analyses 

e Frequency content analysis (FFT) 

REXOR II provides run-time output abilities in printed form and CALCOMP 
plots of selected signals. Deferred signal processing is provided at a 
tradeoff cost, the cost being the peripheral storage required to save the 
signeds for processing later. If the peripheral storage allocated for signal 
data is other than SCRATCH (i.e., disappears after the run), then an economic 
recovery mechanism is automatically provided if the central computer should 
fail during execution. 

It should be noted that the concept of deferred processing does not have 
to be restricted to the type of signal analysis. The processing environment 
is also optional. A time-sharing computing environment with foreground access 
to the saved signal set presents an added dimension of effective use of infor- 
mation to a costly c<miputer model such as REXOR II. 
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2. MAIN - LEVEL 1 EXECUTIVE 


The MAIN program is the REXOR II level 1 executive control program. The 
level 1 hierarchical chart presoited in Figure 2 shovs the separation and 
collection into input-process-output functions. The input routines include 
actual data-read routines, as veil as certain data preparation and initializa- 
tion routines. A flow diagram is presented in Figure 3. The loop structure 
depicted provides multiple case processing. The determination of a case to 
be processed is made in the input routine, READ IN. The software development 
presrated here is in relation to an IBM Installation. Slight variations are 
required to effect the CDC overlay structic'e. These are discuesed in 
section 2.7. The generalized coordinate systm is defined for each case by a 
small set of input indicators. These indicators and required relationships 
for program control are described below. 


Generalized Coordinate Selection 

The eloaents o. .ue generalized coord vector are determined by the iser 
via input. The total syston consists of subsystems which will be identified 
with some physical entity. 



Each subsystem can be optionally included by projjer identificati-'n at the sub- 
system level or in some cases, at the coordinate level, i.e., the selection of 
coordinates within a subsystem can be made. Required inputs and associated 
relationships are presented below. 

2.1.1 Blade subsystem 


• Input 

NB N, Nvnnber of main rotor blades 0 s H, s 7 

D ° 

=1 indicates presence of bending mode #1 
=0 indicates absence 

IM2 Bending mode #2 indicator 

IM3 Bending mode #3 indicator 




Figure 2, - REXOR II - level 1 hierarchical chart 
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NBP Pitch horn bending indicator (elastic feathering) 

NPT N(t Dynamic torsion indicator 

T 

• Relationships 

Total number of bending modes per blade, 

NMB = IMl + IM2 + IM3; 0 < NMB < 3 
Total number of modes per blade, 

NM = NMB + NBP + *IPT 

0 ^ NBP + NPT £ 1 (Mutually Exclusive) 

0 £ NM £ h 

Total number of coord in the blade subsystem, Ngs 
NBS = NB • NM 0 £ Ngg ^ 28 

2.1.2 Svashplate subsystem 

• Input 

NSP Ngp Number of svashplate coord. 

• Relationships 

Specific values identify specific subsets. Trivially, 

Nsp = 0 No svashplate generalized coord. 

Continuing, 



IT 



P.1.3 Reference subsystem 


• Input 


NREF N, 


Number of REF coordinates 


{w} - {*• *• 


2.1. U Sh:.ft subsystem 

All 6 shaft coord, |x, Y, Z, 4', 0, are independently selectable. 


• Input 


X„ indicator 


A v£ilue of 0 means exclude 


A value of 1 means include 


• Relationships 

Total number of shaft coord, N. 


NS = ISl + IS2 + IS3 + ISU + IS5 + IS6 


0 < N„ S 6 

D 


2.1.5 Rotor subsystem 
• Input 


Number of rotor coord 0 £ N^ s 1 
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Total number of generalized coord, 


mOTL 


NBS + NSP + NS + NRLF + NR 


0 S N 


TOTL 


< kk 


Subsystem vector locations within generalized vector are; 


LBS 

^S 

1 

LSP 

^SP 

Ss * “^BS 

LS 

T — 

^SP * ”sp 

LREF 


-S ^ S 

LR 

= 

^EF * ^REF 


A summary of required inputs and the associated relative addiess (RA) 
nvaibers is given below 


Input 

RA 

NB 

1 

I^a 

2 

IM2 

3 

IM3 

U 

NBP 

5 

NPT 

6 

NSP 

7 

NREF 

8 

ISl 

9 

Io£.' 

10 

IS3 

11 

ISh 

12 

IS5 

13 

IS6 

lU 

NR 

15 
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2.2 Input and Initialisaticn 

2.2.1 HkTAF . > Subprogran DATAP is a data>preparatiaa routine. It perforas 
three functions. 

• Lists the eoaplete data deck 

a Prepares a data strean for the input alisoritla of subnrograa REABIS 

• S«ses and processes a special data block. 

Hie first itea is self-eirolanatory* Tbe input data deck is listed as input 
except for the special data block to be discussed. Ibe second function, to 
{urepare a data streaa, requires further cotn t. The input routine READIX 
reprocesses an input record sore than once. Xechaniaas to perfora that task 
are traditionally installation and softvare deponent. DATAP builds a teta 
stresB with nultiple cc^ies of a gim record where needed. 

A xuirsMl data deck my be p .ceded by a special data block containing 
aain rotor aerodsmasde data. This data is handled bf its own set of routines 
iddch are discussed in Secti<»i 2.3. 

A flow diagraa of DATA? is presented in Figure k. It should be noted 
that the input streaa is the standard FORTRAX input laiit 5, and the pr^Mred 
data stress is FG^TRKl unit 12. 

The subprogrsB azgiaaent, XUN^T, an inforaational paraaeter which 
passes inforaatimi fron the aerodynaaic data input routine ABLKX back to 
NAIX for use in other aero routines. This is discussed further in Recti'^ 2.3. 

2.2.2 READIX . - Subprograa READIX controls the data deck processing as 
described in Section 3.1. of Voluae III. This processing includes data unit 
recogniti<xi; case data construct!^; the printing of naster data and case 
data; the punching of an edited naster data unit; setting signal, plot, and 
print default ealues; ai^ finally the processing of signal label and title 
cards. 


Argunents List 

XDATE 8-character date, input 

CASEID 6-character case identificati<», output 

SAVED Logical variable for control 

.F. - naster data not saved 
.T. - naster data saved 
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5CMDRE 


Logical variable for control 


.F. - nore data to process 
.T. - no oore data to process 

CASh LogicEil variable for control 

.F. - not a case data unit 
.7. - a case data unit 

A flow diagraa is presented in Figure 

2.2.3 DTPREP . - Subprogras D7PREP is a data preparatitHi routine. Hain rotor 
blade bending node data is asseobled for the aodes desired. Mode data for up 
to three aodes &ay be defined in the input addresses. Ihe aodes desired can 
be any sul^et of the defined aodes. Also, aodal data is defined in a left 
handed Z down coordinate systea. This routine rotates the data into a Z up 
system. 

Data associated with the shaft subs/stes aay require shuffling also, if a 
coordinate sul»yst« is chosen. Shaft spring and da^wr input aatrices are 
rev and coluzm shuffled accordingly. 

2. 2. It IMIT. - Subprogras IHT perforss a collection of loosely bound func- 
tions associated with the initialization of a case. These functions include 
the one time calculations of subsystoi inertias and oass aatrices. The units 
of sosae ini»its are converted. Certain transformation matrices td:ich will 
main cemstant tiirou^out the case are calctdated. Finally, the generalized 
coordinates and the auxiliary coordinates are initialized. Certain equations 
of interest are presented below. 

Main rotor blade mass: 


Blade inertia: 


”bl 



c(S}ds 





I (S) ds 

^LE 


a(S) (S) ds 
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Fign 


2k 



3- - Concluded 







Total inertia: 


Naas of the transaission mount relative to the huh 

N. ■ N'N 

where [MoJ is given hy equation 351» Volume I and [Ho>p3 siveo equation 
333* Volume I. Also, the fuselage mass [Mo3pt given hy equaticoi 362 of 
Voliaw I, is calculated. Transforms calculated in IHIT include the transform 
from rotor to blade. For each blade, n, 

fcos sin 0 

sin *^oos 

0 0 1 

where 

♦bl ' (n-l)/»^. 

“"n 

The fuselage-to-transmission mount base is also computed 

[^F-Mb] " ^♦Og* ®0g* ^Og^* 

In the notation above, Tp, depicts the standard Euler transformation. 

The remaining initialization details will not be repeated here. The 
reader should refer directly to the code. 

2.3 Aerodynamic Data Input and Preparatioo Modules 

The main rotor aerodynamics of REXOR II are handled via a table lookiq> 
method for the functions: 

S * S 
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A representative set of data is built into the program. Bovever, the user 
may provide his ovn tables to the program. Puzi;her» two sets of tables may 
be provided vith a provision to shift from one set to the other as a function 
of blade station. Another requirement of this simulation is the compatibility 
of input format with those of the C-8l Rotor Program. Compatibility has been 
achieved by installing the necessary input and interpolation routines borrowed 
from the C~6l program. C-8l routines installed intact or in a slightly modified 
form include: 

ABLOCK 

REDATB 

TAK)UT 

TABFIX 

CLCDCM. 

The function of these routines and others will be briefly described below. 

2.3.1 ABLOCK . - The user can supply to the program on a <aice-per-run basis, 
either a replacement set 2 of aero data or an added set 1 to the existing set 
or both. The details of this mechanism are described in Volume III. If an 
aerodynamic data set is sensed by the routine DATAP, then ABLOCK is called. 

This routine controls the input of the aerodynamic block sensed by causing 
the data to be read into the proper tables via subprogram RHDATB. 

Arguments List 

IN FORTRAN input \mit. Designated as 5. 

SETID Set I.D. as sensed cud supplied by DATA?. It is an 

alphaniuneric string of the form 

' 1 • 

A A A 

or 

• . 9 * 

A c /V A 

NUMSET An identifier which will indicate which sets have been input. 

=1 set one only 
=2 set two only 
=3 if both sets 

2.3.2 REDATB. •> Subprogram REDATB simply reads in an aerodynamic table. 

All tables, i.e., C^, Cq, and Cn are of the same form. Subprogram ABLOCK, 
which calls REDATB, controls which table is read. The table formats are 
detailed in Volume III. 
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Arguaents List 


C 


nc 

vz 


1 


Table oaae 


Table size paraawters 


IM FDRTRAH input unit 

2,3.3 TAB(XJT. - ^bprograa TABOUT performs two functions. First, if tables 
have been input, TABOlfT will list those tables on the standard output device. 
Second, all tables to be used by the progran are further prepared for the 
interpolation roxitine. This preparation is initiated for the required 
tables by calls to routine TABFIX. 


2.3.^ TABFIX. - Subprogram TABFIX scans the table specified and prepares 
associated information tables for later use by the interpolating routine, 
CLCDCM. The intent is to provide information which will speed up the argu> 
ment bracketing process. 


2.3.3 TABGEB. - Subprogram TABGEH uses the interpolation routine, CLCDCM to 
prepare two functions for use by the dynamic STALL routine. The two functions 
are 


<=L ' =L<“-">U=C 


3C- 

it - 1 a 


=0 


The functions will be in table form suitable for fast interpolation by 
routine XTRPl. 


2.U Output Modules 

2.U.1 PRIHTD. PRIHTC. - This multiple entry routine produces a formatted 
listing of the program inputs. The entry PRIHTD produces a listing of master 
data and PRIHTC produces a listing of case data in the form of an exception 
report. 

2.U.2 OUTPUT . - Ibis routine prints block time histories. The block size 
is up to six rows of eight variables each. The quantities to be printed are 
selected from the signal set with the print pointers, PP(U8). The corres- 
ponding titles are gathered from the PLAB array, and the first eight charac- 
ters of each signal title are used. The first and last, points are always 
printed. An input, PFREQ, controls the print frequency. The user control 
of the block time history output generated by this routine is discussed 
in Volume III. 
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Arguments List 


SIGUNI The number of the FORTRAK unit which holds the signal data. 

The value is set in the calling routine, MAIN. 

NPTS The number of points saved in the signal set. 

CASEID 8-character case identification. 


2.U.3 HARM. - This subprogram will provide harmonic analysis output for up 
to 50 signals. The output will consist of the Fourier coefficients 


a 

o 


a 

n 


b 

n 



f(x) dx 


f(x) cos nx dx 


f(x) sin nx dx 


for n = 1, , 8 (8 harmonics). 

The data will also be presented as phase-amplitude information where: 



If b < 0 
n 

= 0 + 360/n 
n n 

Subroutine limitations are arbitrary and include: 

• The period of the functions analyzed is 2ir 

• Eight and only eight harmonics 

• A maximum of 50 functions can be analyzed. The functions must be 
defined in the signal set. 
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Arguments List 

NPTS Number of points in signal set 

SIGUNI The number of the FORTRAN unit iriiich holds the signal data. 

CASEID 6-cbaracter case identification fc* output 

TCYCLE Numbs: of cycles consumed in TRIM 

2.1*.U FORIT. - Subprogram FORIT is derived from the IBM Scientific Subroutine 
Package. It performs a Fourier analysis of a periodically tabulated function 
by coB^uting the coefficients of the desired number of terms in the Fourier series 
series (Reference 1): 

M 

F(x) ~ *o * ^ \ ^ ^k 

k=l 

Arguments List 

FNl Vector of tabulated function of length N 

N The number of samples over the interval (O, 2it) 

M Maximum order of harmonics to be fitted 

A Resultant vector of Fourier cosine coefficients of 

length M+1 

B Resultant vector of Fourier sine coefficients of length M+1 

lER Resultant error code vrtien 

lER =0 no error 

IER=1 N<2M+1 

lER = 2 M < 0 
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2.U.5 GRAPH . - Subprogram GRAPH controls the generation of CALCOMP plots. 

This includes the opening and closing of plots as required by CALCOMP software. 
Further, a logical decision is made as to whether there are enough points to 
plot. Plotting is not accomplished unless sufficient points are available. 

The plotting of signals is performed by a subordinate routine, RCPLOT. 

Arguments List 


NPTS Number of points in signal set 

SIGUNI Logical unit number of signal set 

MORE Sign-off signal. A logical variable idiich if fedse then 

sign-off. 

FENTRY First entry signal. If .TRUE, then first entry, so open 

plots . 

CASEID 8-character case identification. 

A diagram is presented in Figure 6. 

2.k,6 RCPLOT. - Subprogram RCPLOT is the CALCOMP plot generator. Up to 
100 signals from the signal set can be selected for plotting. Only a few 
words concerning plot specifications are presented here. Details concerning 
user options and control are set forth in Volume III. 

Although the nximber of points saved in the signal set is open-ended, 
dimension statements within this routine have arbitrarily set the maximum 
nimber of points per signal at 20 't 8. This nmber can be extended. 

Plotting proceeds in a frame-oriented manner. Signals are collected 
four at a time for plotting on a frame. A minimum frame size corresponding 
to a standard 8-1/2” x 11” page is produced so that CALCOMP continuous 
paper rolls can be cut page size. If time historji' signals ere longer than 
the minimum, then the abscissa is continued the necessary amount. A sample 
frame can be found in Volme III. 

Argiments List 


NPTS 

Number of points in signal set 

TSCLE 

Abscissa scale factor, sec/in. 

NFREQ 

Plotting frequency, i.e.. 


=1, plot every point 
=2, plot every other point 
=3, plot every third point 
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SIGUNI 


Logical unit containing signal set 

CASEID 8-character case identification. 

2.1*. 7 BSCALE . - Given a signal to be plotted, BSCALE determines a plot scale 
factor such that the signed, is restricted to a 2-inch channel. 

Arguments List 


X Array to be scaled 

N Number of points in array 

DY Scale factor of the form 

DY = 1 X 10^*® 

DY = 2 X lo’^^ 

DY = 1* X 10^^ 

NA Exponent above . 

2.1*. 8 PDATE. - PDATE is an installation-dependent routine designed to provide 
the date in cheuracter form. The subroutine argument is an eight character 
variable of the form: 


NDATE = MM-DD-YY 

2.1*. 9 PAGE. - Subprogram PAGE performs a page eject on unit 6 idien called. 
A page header is produced with the case identification, date, and REXOR II 
input titles 1 and 2. On exit, a line coxinter is set to 5. 

Arguments List 


ICT Line count on exit 

CASEID 8-character case identification 

NDATE 8-character date 

TITLES REXOR II title array with dimension (15, 1») 
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2.5 MASGEN 


Subprogram MASGEN is an alternate process control program. It is 
optionally invoked from MAIN. Its function is to generate a generalized 
mass matrix via perturbation techniques where 



''■'k ' - “’’k 



A mass matrix generated by this method provides an excellent diagnostic tool 
for studying the correctness of the computed mass matrix ^en new equations 
*ure developed. MASGEN is not necessary to the operation, but is provided for 
possible fut’ore use. A flow diagram is presented in Figure 7. 


2.6 Block Data Subprograms 

A number of block data subprogrt-ms are present in the program. A list 
of the subprograms and the LABELED COMMON'S witiiin each is presented. 


NAME 

COMMON LABEL 

DESCRIPTION 

BLOCKA 

/TAB / 
/CLACLO/ 

Main rotor blade aerodynamics 

BLOCKB 

/CONSTS/ 

Constants 

BLOCKC 

/NAMES / 

Input parameter names 

BLOCKD 

/LABELS/ 

Signal labt sfault set 

BLOCKE 

/RAINIT/ 

Default relative addresses for the signal 
set, print, and plot tables 

labeled COMMON /CONSTS/ 

is expanded below. 


NAME 

VALUE 


G 

32.171^ 


PI 

3.11*1593 


TWOPI 

6.283185 


PI02 

1.570796 


DTOR 

57.29578 


RTOD 

0 . 0171*5329 


RHOZ 

0.00237689 


VSZ 

1116 . 1*5 
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2.7 CDC Orerlaj la^loMDtation 


The iapleaentetioB of an ororlay structure in a CDC environaent requires 
the restructuring of the software into a series of prograas. The execution 
flow aost be explicitly defined withiB the source stateaents of the various 
prograas. The changes required to amvert IBN software as described here to 
CDC software is as follows. The reader can refer to Figure 3* REXDR II NMH 
prograa. Subprogram READI3 is elevated to prograa status and placed in 
priaaxy overlay 1. It is invoked by a CALL OTBItLAT (PFILE, 1, 0, RECALL) in 
place of CALL RSADDf . 

The CALL SUSEFl. CALL MASCES, CALL FLT sequence of the MAH prograa are 
reaoved, and given pragma status as PROGRAM PROCESS. They are placed in 
priaary overlay 2. Subroutine HARM is also giveu prograa status, and placed 
in priaary overlay 3. Finally, subroutine RCPLOT is girm prograa status, and 
placed in priaary overlay b. The respective calls have be«i altered. Data 
ci:niiilcatic»i8 idiich existed through argiaent lists are handled through 
labelled CGMKHIS. 
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3. FLY - LEVEL 2 EXECUTIVE 


The IcTel 2 hierarchical diart for FLY is presented in Figure 8. FLY 
is a secondary executire routine which controls the calculation of the 
generalized coordinates, thus, generating all the desired tine signals. FLY 
is prinarily a loop structure built around the integrators as shown in 
Figure 9« The eleaents of the loop are described below. 


3.1 FLY 

A flow diagraa cf FLY is prearated in Figure 10. A nwiber of prograa 
counters and run controls are defined below for further clarification. 

Ihe progran Tariable TRIMED is a logical wariable which is initialized 
as fhlse and under certain conditions is set to the value, true, within the 
TRIM subprogram.. The wariable TRIMED indicates if the Thiele trianed. The 
variable MCCMfP is the generalized aass aatrix cooputation flag. It is 
initialized as true. The variable OISFNX flag is the display natriz flag and 
is initialized as false. Otter variables are: 

KPTS Huaber of points saved in the signal set 

PTIME Trajectory tine seasured froa end of TRIM. Its priaary 

use is as the arguaent to the pilot control tables 

IC I Integration step counters used for logic control within 

ICHTAX I the integration routines. 

AZIMIS Blade aziicuth step counter within (me rotor revoluti<m 

CYCLE A count of coopleted rotor revolutions. Blade 1 is the 

reference 

CTEST Trim convergence test (K>ntrol counter 

TCYCLE A count of the number of cycles required to trim 

LTICK Main rotor pip counter. Counts blades crossing the 

fuselage 

HCC The number of second order differential equations in the 

generalized coordinate set to be integrated. This count 
is sdtered in TRIM 

RE The number of auxiliary differential equations to be 

integrated. 
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One should note that if in the future nev auxiliary differential equati(»is 
are defined, the Tariable. HE, defined above nust be changed to the corre c t 
count, .'t should also be pointed out for future consideration that a fUnda^ 
■ental asfin^ion entry into FLT is that the goieralized coordinate mt 
has been defined and initial ^^ues for the vectors GC(I), GCD(l), and 
GCDOd) have been aade. 

The arguaents in the calling sequence to FI>¥ are priurily output 
oriented. 

CASEID 8 character case identification which originates in 

REAOIH 

lER Error indicator. Initialized in MAIH. Value set within 

ACCEL. 

SIGUHl FORTEAH logical unit for storage of output signals. 

The renaining arguments HPTS and TCTCLE cure defined above. 


3.2 COC«D 

Given, the generalized coordinate vectors, and the necessary indicators, 
this subprogram defines the program problem coordinates. The generalized 
coordinate vector is 
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A0VMKE 

THC.ABMU1M 

AMDCVCU 













siailarly for 

{'} 

and 

1110 foUoving definitions and coord options aire aradlaUe. 
Blade System 




*BS 


For each blade* n 


If • • • t 


Bending Modes 

Nn " Wb 


m « 1, 


• • • » 



Likewise for 
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and 


w. 

Pitch Horn Bending (Elastic Feathering) 
If BBP * 1, then 

Likewise for 


P and P 
n n 


Also, since pitch horn bending and dynamic torsion are mutually 
exclusive. 



0 


Dynamic Torsion 
If NPT = 1, then 


Likeid.se for 


and 


Also 






= p = p 


Ul 


0 



NOTE: Quasi-static pitch horn bending and torsion 

are available as adternative models. 

NOTE : No provision has been made at this time 

for the special case where N6S ^ 0, i.e., no 
blade system. 

Svashplate Options Ngp = 0, 1, 2, 3 
Ngp = 3 All active 


*SP = *^1 SP 


®SP = *^2 5? 
^SP = ^3 SP 


Similarly for 


♦sp *SP 


6cp 9; 


SP 


^SP ^ ^SP 


“sp’ " 


(^gp and 0gp as defined in Option (®gp ® 3) 


^SP * ^SP " ^SP ^ ° 


"SP = ^ 


A * A ** 

^SP ' *^1 SP’ ^SP “ *^1 SP* ^SP “ '^1 SP 

ri= 


fcos 4) -sin 4» 

' c c 


9SP 

/^SP 


i l’SPl 

®Spl 


sin cos 41 

c c 


cos 4/ -sin il' 
c c 

sin 4< cos it 
c c 



(-■See *c) 

< 

! >^0 

- 

rtc *c ] 

i 



h2 



0 



*Scc “ *c 

BA (123) 

*SfC “ 

RA (12U) 



Zgp is defined as in option (Hgp * 2) 

(^gp and Ogp are defined as in option (Ngp 
Reference (l*u£p “ 6) 



Likewise for 




and 


1 ) 


»*3 




REF 


Shaft Subsystem Options Ng = 0, 1, 2, 3, *», 5, 6 



Likewi se for 

{^} 

and 

N 

For options where Ns < 6, a special vector compress-decompress routine 
assigns the appropriate coordinates. 

Rotor Subsystem 

^R * {‘Ir} 

= {^r} 

^R {^r} 


uu 



0 


'*'r “ 

ii»_ = INPUT 
n R 

♦r ■ ^ 


3. 3 CMPRES 

This siibprogram will compress a vector, X, to a vector, Q; 

{*} H 

(6x1) (Nxl) 

or decompress a vector, Q, restoring X, 

{"} 

(Nxl) (6x1) 

by using t mapping function, IND, of indicators. The indicators operate 
positionally. If an element is 1, the corresponding element of X is passed 
to Q. Also, 
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i=l 

Descriptio. >f Arguments 


Input 

COMPR logical variable 

If COMPR = .TRUE., then X is compressed to Q. 
If COMPR = .FALSE., then the vector X is 
reconstructed from Q. 

IND Vector of indicators used in the mapping. 

INDi *» 1, then element i is to be passed 
INDj * 0, then element i is to be ignored 
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Input /Output 


If COMPR = .TRUE. 

then X(6) is input 
Q{N) is output. 

If COMPR = .FALSE. 

then Q(N) is input 
X(6) is output. 

Undefined X's will be set to zero. 

3.1* TRIM 

Subprogram TRIM provides a means by which the equations representing th 
vehicular system can be put in a so called state of equilibrium. Actually, 
a more general capability is available. Find x such that 

g(x) = f(x) - c = 0. 

The algorithm to solve this problem is simply stated 

x*^ - k • h*^(x) • dt 
h’^'^Cx) + k (f’^(x) - c). 


n+1 


h^U) = 


The user can choose up to ten functions and ten corresponding independent 
variables from a supplied table. The collection of functions is not treated 
as a system but rather as independent relationships. Convergence is deter- 
mined not by monitoring the function, g(x), but rather the relative change 
in the independent variable over one rotor revolution. Convergence is 
assumed if 


X - X* < e 

for all X. X* is the value at the end of the previous cycle. 

A general, not detailed, flow diagreun is presented in Figure 11. 

It should be noted that during trim the reference subsystem and the 
rctor subsystem are not integrated. Thus, as indicated on the diagram 
certain coordinate overrides must be performed. 




TRIM BACKUP HAS 
BEEN REACHED. 
SIQNAI TRIMMED 



rPIM flow diagram. 
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The coordinate overrides are 


# = input 
£» 

®E “ ^ ^ 



fg = f(Lp), Lp = input 



♦p = Ot 


Other inputs idiich are central to the tria process are listed belov. 


Inputs 

Description 

HA 

BC 

Humber of contn>l lavs 
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CVP(IO) 

Control variable pointers 

1301-1310 

FUliCP(lO) 

Control function pointers 

1311-1320 

CVE(iO) 

Convergence errors 

1321-1330 

CVK(IO) 

Control variable gains 

1331-13UO 

FUHCD(IO) 

Desired function values 

671-600 

TCUI 

Tri* cutoff, cycles 
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Tbe description of arguaents are: 

fSIMED Logical variable indicating trim condition 

CYCLE Main rotor cycle count 

CTEST Tria convergence cycle indicator 

IC Generalized coordinate integrator control 

’iCC Generalized coordinate count 

TCTCI£ Cycle nuaiber vfaen s y s te m declared trijawd 


3.5 RETRVE, ST(»{E 

This is a Bultiple entry subprograa designed to retrieve and store trin 
paraaeters viiich are defined in various COMOIS. The process is table 
driven. The current tria paraaeter definitions are givaa in the foUoving 
table. 













Parsuseter 

K.D. 

HuBiber 


Ccmimon Block Ident 


Block 

Offset 

Blank. 







ll» 


n 

X 





63 

**ZZejid 

15 


X 





92 

’b 

16 

X 






320 

**ZR 

17 

X 






963 


18 

X 






961 

**YR 

19 

X 






962 

^ZR 

20 

X 






960 




The retrieve/store aechanisB is depicted in Figure 12. 

A description of the arguments is 

H Hiober of paranieters to retrie/e/stoire 

POUTTR Vector of H pointers. The values are drawn from the ID 

nxaabers defined in the definition table. 

V Vector of parameters retrieved or stored. 


3-6 TFMIM 

This program cooputes all tine varying transformation matrices. 
Input Data 

All necessary singles are in blank coonon 

Quantity Haae 

PSIR 

PHSP 

THSP 

^ PHIS - DISPS(U) 

9^ THTS = DISPS{5) 

PSIS E DISPS{6) 

# PHIE 5 AUXE(l) 

9!: THETE 5 AUXE(2) 

PSIE = AUXED) 
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N.^ NB (XMtON/DATASl/ number of sain rotor blades 

Tp TRBL(3,3,7) CO»#ION/TFOItt6/ transforas from rotor to blades, 

n coB^uted in MAIH. 

Subprogram TFORti uses one subordinate routine, EULEBT. EULEBT evaluates 
a general Euler transformation. The ose of EULERT is expressed functionally 

as 


Tg (♦, e , f ) 


Output 

All transforms appear and are available thru labeled (XMOH/TFOBMS/ . 
The transforms are defined below. 

Hub to Rotor 


r 1 

-cos fg 

slnTp 

0“ 

r - “J ” 

Sin Tg 

CO. 

0 


0 

0 

-1 


Hub to Blade n 


RbL — H " "^BL — R "^^R — H 
I n n 


for each blade n = 1, ...» Ilj, 
Hub to Svashplate (Rotating) 



* ®SP* 


Mount Base to Hub (Transmission Isolation System) 


jj^S ^ hJ 2 I^MB — hJ * ^E 
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Reference to Earth 


['*^REF — eJ 


•^E ^-♦e* -®E» - *E^ 


3.7 EULBRT 

This program calculates the Euler transforaation matrix given by 


where 


Tg « T(^) • T(e) . !(♦) 


T(*) 


cos<> 

-sint 

0 


sin'l' 

cos’i' 

0 


0 

0 

1 


T(0) 


T(*' 


cose 

0 

-sine 

0 

1 

0 

_ sine 

0 

cose 

1 

0 

0 

0 

cos ^ 

sin# 

0 

-sin ♦ 

cos# 


Input 

« - PHI , Radians 
6> THETA, Radians 
’I'- PSI , Radians 
Output 

T- - TE(3,3) 3x 3 transformation matrix 

£ 
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3.8 MOTIOH 


Given the baisic reference TOtions, this routine computes all najor linear 
and angular motions except hlade and blade elesMnt w>tions. They are computed 
elseiriiere. The equations of interest are presented below. 



Referoace Acceleration 







Earth Motions 



and 


AUXD(21) = Ug 
AUXD(22) = -v_ 
AUXD(23) = -Vg 


Fuselage 




P 

4 

f 


REF 






FORTRAN 
* VELRF (4) 

= VELRF (5) 
= VELRF (6) 

= ACCRF (4) 
= ACCRF (5) 
= ACCRF (6) 

= AUXE (4) 
= AUXE (5) 
= AUXE (6) 


55 



= ACCRF (1) 





* ACCBF (2) 
= ACCBF (3) 


Velocity and acceleration of MB relative to fuselage 



use routine LINVEC vith arguments 
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Mount Baae (Transoilssion support) 



VDMB^ 




1 

' ^ 1 

1 ( 

p \ 


4 

: ' K 1 ' 

I 4 

1 



\ 

1 

^ 1 

Imb I 

r 


F 
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Hub 


Compute a set of intermediate vectors 


|ViX:[ and {VEc[ using LINVEC with 





X 

Y 

Z 


'MB-H 


X 

z 


.REL 


fMB-H 


Y 

Z 


kREL 


MB-H 
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Using program AKGVEC* compute 



The arguments are 
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Rotor 


Quantity FORTRAN 




Svashplate 

Quant ity 




= WR(1) 

= WR(2) 

= WR(3) 

= WDR(l) 
= WDR(2) 
= WDR(3) 

FORTRAN 
= WSP(l) 

= WSP(2) 

= WSP(3) 


i pj = WDSP(l) 

4 / = WDSP(2) 

r) gp = WDSP(3) 


VDR(l) 

VDR(2) 

VDR(3) 


Eq 7h of Vol. I 


Eq T5 of Vol. I 


Eq 159 of Vol. I 


Eq 160 of Vol. I 
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Gravity in the rotor systi 


{*} R “ h - h] h ♦ mb] [’iffl f] {*}ref 

3.9 LUVEC 

This s^a>p^osra■ ctmp\xtes general linear inertial velocity and accelera- 
tion vectors of the font 

{»} . {v^ . . [.] {vj 

- {'4 * {’1} * [“] H {’1} * M {»i} * " [»] {^} 

where 

0 -W(i) ¥(2) 

« y(3) 0 -W(l) 

-V(2) W(l) 0 

and siBilarly 


c 

-H(3) 

W{2) 

[v] * W(3) 

0 

-w(i) 

-W(2) 

w(l) 

0 


The subrou^-ine arguaent list is 

(Vi* ir^. Vg. Vg. w, v, v. v) 
all arguaents are (3x1) vectors. 
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3.10 AIGVEC 


The subroutine computes goieral angular Telocity and acceleration vectors. 
The arginsent list is (HA. HDA. A. AC. ADD. VB. VDB) vfaere 

Inpits 




defines transforaation froa systea 
a -*• b 


derivative of A 


% 

derivative of A 
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Outpjt 



Algoritha 




0 

0 

T(#) = 

0 

cos # 

sin # 


.0 

-sini 

cos^ . 


"0 

0 

0 

T(*) =$ 

0 

-sin $ 

cos # 


.0 

-cos • 

-sin ♦. 


0 

O 

6 0 

« 

-sin d 

T(e) = 

0 

1 

0 


jsin 

e 0 

cos 9 - 


6U 



0 -cos 6 


T(e) » 


T(») « 


T(<^) - 


* 0 



* c 

ly 

\ c 


'-sin 8 
0 

cos 8 

cos <1' 
-sin <t> 
0 

-sin t 
-cos ♦ 

0 



0 0 I 

0 -sin 8 J 

sin ♦ O' 

cos 0 

0 1 

cos 0 

-sin f 0 

0 o: 
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Finally, 


and 




°\ • 

0 > ♦ T 


“I 


T(ii-) {vg} ♦ T(*) I Vg| 


3.11 CHTROL 

Subprogras CBTROL assembles pilot controls with trim inputs and stability 
augmentation iiq>uts. Ibese cabined inputs then operate the rotor and fixed 
surface panels. Controls are presented in Section 7 of Volume I. The general 
flov of tte routine is depicted in Figure 13. Kotice in the diagram that an 
abbreviated set of equations is used dxiring TRIM. Whereas, during FLY, pilot 
control tables and SAS inputs are active. 

During trim, the control equations are givu by the foUoving equations. 


^o ■ ^eFS 


variable 

^p ■ *S>FS * 


variable 

DB ' *4)BFS 

DB,IC^ 

constant 

FL “ ’SfXFS 

^ FL,IC^ 

constant 

iw iwFS 

iw,IC^ 

constant 

HT “ *^iHTFS 

^ iHT,IC^ 

constant 
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6 

6 

0 


'el 

AIL 

RUD 

°TR 

X 


Y 

c 


* *^AILFS 

* *SlUDFS 

“ *Sprfs ‘ 

" *ScCFS 
~ *SfCFS 


°EL * *c.T^ ■*■ 'eL,0 


^AIL ■ ^c,T^ ■*■ "^AIL.O 


SuD * *‘c,T^ * *RUD,0 

(i * (f r _) + S ♦ 

^ °ra,o 


'"c ' 

® 300 


®o 

300 


^^OTR 

36o 


»o 


During FLY, a number of pilot controls may be input as a function cf time. 
The controls available are Xcp, Ycp, rcp, 5pp, 8op, 5QBp, 5FLp, 6iVp, 5iTOfp. 
These functions are evaluated via linear interpolation. Continuing, SAS inputs 
are available for six functions, Xcgj^, Ycgj^, SEIsaS* AALgAS* ®^^SAS* 

SRU%AS« ^AS logic is found in routine, SAS. The control equations are 

straight forvard and can be found in Volume I and in the program listing. 
Finally, stick actuators are applied to the longitudinal and lateral stick 
controls. Actuator logic is depicted in Figure lb for the X-stlck. Similar 
logic applies to the Y-stick. 
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‘ic ” irPiM 


IWOCHANNELS X. 

AND 

Yc 

iNrUTS 

MAKE 

FORTRAN 

RA 

^ M 

TAUXC 

292 


TAUYC 

293 

f MVC 

RATLXC 

437 

*i(rrc 

RATLYC 

438 


Figlire lU. Actuator logic. 
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3.12 SAS 


This is a general SAS control routine. The SAS block diagram is pre- 
sented in Figure 15. The digital implonentation of that diagram is 


Si . . Kj • ♦ II 




S3 - Si - 


Su = "1" <|S3l- • Vl^sl 


I, - t • S; 


h‘h- h 


K. « — 

** ^1 


^SAS “ ^ 


All integrals have zero initial conditions 




= 0 


The error signals are of the form 


^ ■ q - q^RiM 


The routine is invoked by the call 


CALL SAS(£A, EB, SASD, II, 12, IDl, ID2, XSAS) 


where 


EA 

EB 


a 



I 

I 


Error signals 
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(ANALOG IMPLEMENTATION) 


Figure 15. - PAS block diagram. 






SASD 


Vector of constw.ts for the particular control where 



XSAS 


Feedback gains 

Filter gain 
Frequency breakpoint 
Output gain 
Filter time constant 
Limiter 

Low pass filter integral 
Washout integral 

Integrands to be integrated elsewhere 
Control output 


3.13 CYCLIC 

Subroutine CYCLIC provides the linkage between the swashplate and the 
main rotor. 

The main rotor collective is determined in CNTROL, However, the asso- 
ciated derivatives are determined here. 


e 

o 


= - 


0 

o 


= - W' 


The cyclic angles are 





1 


T2 



IS 


• (i) '*8? *0 * «SP 




®1S ' ItJ ‘*SP ■=“= ®G ■ ®SP “o' 


The derivatives are 


*18 “ (e) '♦sp ‘a * »sp «o' * (f) j «» 

® 1 S ' (i) ’♦sp ®G - «SP ®g’ * (l) J 6 o 


*18 - (e) 


[Vgp Sin Sq + 'e'sp cos Bg] 


^2 (f-j ^ Ugp ♦ Ogp cos J 




+ lei 0 A-„/(d/E) 

O lb 


B 


: , ± 


Is ■ Ve y ^*^SP ®G ■ ®SP 


+ 2(-f-| L cos 3„ - 6ep sin S, 


6 o '<^SP '^G ■ «SP ^^G 




+ lei 0 B_/{d/E) 

o IS 


A^g/(d/E) 

B^g/{d/E) 
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3.ll* PDEPIV 


This routine computes partial derivative matrices relating a set of 
dependent coordinates « to a set of independent coordinates, tj. Not all 
program partials are computed here. Those computed in this routine are 



where 



7U 



3ts 



• 

-►hJ 



cos 4'. ' 0 sin 0 

s s s 

0 

-sin t cos cos 0 

5 S S 

- 

sin e 0 1 _ 



•^5 




[*^F mb} : 

JL.^ 


(’'ref) 

[ “ 

• 

• 

• 

• 

• 

pF - MB 



' 0 

2 -Y 

■ [ 

*f-mb] 

-2 

0 X 



. T 

-X G 

IlfL) 

II 

\T 

aj 


J ^^REF j 

[ 1 

I^^REF 1 


3-15 AUXOOT 



This routine perforss auxiliary calculati<sis for oxitput. The quantities 
coBpoted are stored iu the vector, OUTD, idiich is defined in OOMOlf and is 
therefore availahle for si^al definition. Definitions of the currently defined 
eleaents of OITTI) follow. 

Speed in the reference system 


OUTTd) = ♦ w^„„) 

\ ref ref ref/ 


1/2 


= V_ 


REF 


Angle of attack 


OUTD(:) . su-1 [•'EEr/f ' 


1/2 


( REF * ^ REF ) 


]• 


57.29578 


Angle of sideslip) 


OUTDO) = sin’^ ["^REf/^j^] * 57.29578 


76 



Flight path angle 


OUTD(U) 

OTm)(5) 

0UTD(6) 

0«TD(7) 

OUTD(8) 

OUTD(9) 

outdCio) 

OUTD(ll) 


sin 


-1 




REF 


= AXS * 57.29578 


= • 57.29578 




57.29578 


Attitude gain 

0UTD(12) = h - h (input) 
e o 


Seven functions at three blade stations on blade 1 a:e made available 
prjjuarily for hamonic analysis. The routine will assemble data at the 
stations closest to 25, 50, and 75 percent blade radius. 

The seven functions are 

e Span force 

e Inplaue shear 
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• Flap shear 

• Torsion 

• Flap ■cment 

• Inplane aoaent 

• Torsion, aero only 

I%e functions are equiTsl«iced to the QUID vector, 0UTD(13) through 0UTD(33). 

3.16 SIOCZI 

SIGGES is the signal generator routine. It vill collect, up to a M a xigu n 
of 200, signals froB blank CCMOH. This is accosQaished by a pointer array, 
SLOG, which is input. &ch signal collected can be scaled with a unit cchi- 
verter, one for each signal input. These oultipliers are also input, UC. 

Arguaents List 

STIME Trajectory tiae 

KPTS SiBiber of turn points output on exit 

SIGUNI Logical unit on which the dUita is written 

3.17 DHASH 

Subprograa DtfASH calciilates main rotor and tail rotor inflow functions. 
3.17.1 Main rotor equations 






1/2 



3 


2P(B»R)^ 


- vw_ 


■•MR 
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Pt - 
Sir ■ 32 B B 


k ~ ?T 

p(B-R) Tffl 



>5 

32 B 



(B »)'' 


The derivatives are integrated locally using a Euler integrator 




♦ 



dt 




dt 




♦ w. 


Sir 


dt 


3.17.2 Tail rotor equations. - The tail rotor functions are hrndled hy first 
order log equations. 

Inflow 


n+1 

TR 


= w. 


-dt/x 


+ w. 


TR 






^~dt/ T j 


Longitudined Flap Angle 


TR 




The superscript 
Wi and Ai 

TRjj TRjj 


n, stands for iteration nunber. The 
are coarruted in subprogram LOADS. 


tai.. rotor functions 
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3.18 TORS 


Quasi.-static torsion is optionally cosiputed ly subprogran TORS. The func- 
tion evaluated here is described in detail in Section 5.6.5 of Voluae I. 


3.19 PHORH 

Quasi-static pitch bom bending is optionally conputed by subprogran 
PSORX. The feathering angle incresent due to pitch bom bending is determined 
by a first order log equation for each blade 


idiere 






The reader can refer to Section 5.6.6 of Volume I. 


3.20 ETORQ 


An engine can be modeled as described in Section 5.10 of Volume I. ilie 
engine equation- are evaluated within subprogram ETORQ with the primary output 
being . The pertinent equations are 




. a _ . n-1 
■^GEN ■ '^GEN 




• n-1 


^2 ^'‘'eNG “ %G,TRIM^ ■ '^GEH 


dt 


«XA 


^”eng • 


ENG 


^NG, 


TRIM 


3%EN 


^”erg . - . . 

^^ENG ■ W.TRIm' 
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3.21 IMTGGC 


This subprograa prorides the integration of the generalized coordinates. 
The integrator is an Adaas-Bashforth Wpoint open relation. An accunulation 
of back value inforaation is required. Therefore* a Euler starter is used. 

It can be seen in Figtire l6 that accelerations and velocities are separately 
integrated. This is for bookkeeping purposes. 

A description of calling sequence argwents follow. 

Hk Hiod>er of integrators to use 

DT Integration interval, sec 

IC Step counter. It will be increased by one on ret'urn. 

V Displaceaent vector. Will contain iq>dated values on return. 

Y. Velocity vector. Will contain iq>dated values on return. 

YuD Acceleration vector. 


3.22 UPTGAX 

Subprograa IHTG^Y provides an integration nechanism for all auxiliary 
integrals. It is identical in design and fom with IHTGGC except that only 
one set of integrals are applied. 


Y 


i 



dt ♦ Y 


i/IC 


3.23 JETTB 

Subprograa JETTB is designed to provide a main rotor blade Jettison 
sequence . 

This subprogram is special purpose in that it is only for a ^bladed sys- 
ten with a specific Jettison sequence. A blade is Jettisoned by removing it 
from the generalized coordinate set, i.e., all generalized coordinates are 
redefined. The generalized coordinate integrator back values are also shifted 
so that integration is not interrupted. The rotor-to-blade transformations 
are edso shifted, and are not redefined. Other blade dependent data affected 
is the quasi-static pitch horn and quasi-static torsion. 
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Figure l6. - INTGGC flov diagram. 


o2 




Operationally, JEITB is activated frcm FLY at the beginning of a pre- 
scribed rotor cycle (JETCYC) neasured fron trim. Upon this first entry into 
JETTB, main rotor blades 1, 3, and b are removed from the generalized coordi- 
nate set as described above. Upon exit, the ntaber of blades (NB) has been 
reset to 2. At a prescribed nvmiber of integration steps later, n = NAZ/3, 
JETTB is reentered and the remaining blades are Jettisoned. At that time 
quasi-static pitch horn and quasi-static torsion vill be turned off. Also, 
all svashplate coordinates vill be turned off. 

The user is reminded that some auxiliary outputs are a function of 
blade 1 data. During a Jettison sequence, blade 1 is redefined. Thus , mixed 
data may be processed. 
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U. ACCEL - LEVEL 3 EXECUTIVE 


The level 3 hierarchical chart for ACCEL is presented in Figure 17. 
ACCEL, as a control program, is simple and straight forward. It directs the 
calculation of the components of the force-nass-acceleration equation and 
evaluates the recursion relation 


..t 

q 



. ..t-At 

+ q 


which is fundamental to developoent of the REXOR II model. Further details 
concerning the subprogram ACCEL are presented in the next section. The func- 
tions GFORCE and hiwSSF depicted in bold-face in Figure 17 are treated later 
as level U executives. GM>.SS is discussed in this section and the matrix 
routines CHOSKT and MXOUT are discussed in Section 7. 


U.l ACCEL 


The subprogram ACCEL computes the generalized acceleration at time, t, 
by applyirg the recursion relation 


..t 

q 


-t-At 

= q 


..t 

Aq 


The incremental acceleration is provided by solving the linear equation 


[m] {ah} . {af} 

The elenents of this equation, the mass matrix and the force vector, are 
provided by subordinate ’•outines. 

The variables in t argument list are: 

CASEID 8 character case title for output purposes 

MCCMP Logical variable which controls the computation of the mass 
matrix 

DISPMX A logical variable used to control the printing of the mass 
matrix 

lER Program error indicator 
The details, .f ACCEL are depicted in Figure l8. 


Qk 




CD 

V/1 


Figure 17 



ACCEL leve] 3 hierarchical chart 






COMPUTE {AF} 


COMPUTE !M] 





Figure l8, - ACCEL flov diagram 




U.2 04ASS 


Subprogram GMASS provides three functions. It determines whether a mass 
matrix is to be conputed and calls MASSSP if required. It shifts the mass 
matrix to working storage. This is necessary because the coefficient matrix 
is destroyed on exit from CHOSKY. Finally, GMASS will list tne mass matrix 
if so directed. All of these functions are depicted in Figure 19. 

The system solution routine, CHOSKY, operates on double precision data 
in the I.B.M. software. Double precision is not required in C.D.C. software 
because of the extended word size. 
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5. a^ORCE - LEVEL U EXECUTIVE 


GFORCE is a fourth level executive control program. Its primary output 
is the generalized force vector, AF, which is defined in COMMON as DELF. The 
elements of GFORCE are depicted in Figure 20. GFORCE and its subordinate 
routines are presented in this chapter. 


5.1 GFORCE 

This subprogram assembles the generalized force vector by optionally’ 
computing the required subsystem fc.*ce vectors. This is best depicted in 
Figure 21. The assembled force vector would look like 


/ 

BS 

NBS X 1 


il''«lsp 

NSP X 1 

= < 

i^^ls 

} N3 X 1 

NTOTL X 1 

1 



f 1 REF ' 

i NREF X 1 


1 t‘S R 

1 NR X 1 


if all subsystems were triggered. 

Subprogram GFORCE heis only one argument, lER. If an error is detected in 
SWEEP, lEE is set to 999 and control returns to ACCEL. 


5.2 Main Rotor Analysis 

Subprogram SWEEP and its initialization entry bWEEPl, computer main rotor 
blade forces and generalized mass elements. To this end ‘ ade elem nt notions, 
partials, and aerodynamic forces are also computed. Thi , s tne m. ^or program 
computation sequence, auid will only be outline! here. References to Volume I 
will be made where required. Other subordinate routines are presented in this 
section also. 


5.2.1 SWEEPl, SWEEP. - The chief output of SWEEP is the blade root forces and 
generalized mass integrals. The forces can be presented as 
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GFORCE Level hierarchical chart 

































vhere the integration is over all "blade stations, k. The above set of 

integrals are generated for each blade n=l, , N„. 

15 

Further, the mass integrals can be presented as 

[«] - 

where the integration is over the blade stations and for each blade. Also, 
the matrix [M] is symmetric. 

The matrices are subsets of the partial derivative matrix defined 
below. For each blade emd station 
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also define the blade forces 



The evaluation of the above partials and forces require the evaluation of 
the blade element motions and aerodynamic forces 
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{ °BLe} , { ®BLE } . { °BLe} . { ^Le} 

Blade eletnent laotions are treated in detail in Section ^.5.5 of Volme I. 
Those equations will not be repeated. ather, program names for certain key 
quantities will be presented. 



Blade element freestrean absolute acceleration 



DDXN 

DDYS 

DDZH 


Blade element angular velocities and accelerations 



ANGBL(3) 



9k 


angbldO) 



Partial derivative matrices are identified as 


where 




PTAUM(3.3) 



TBLR(3,3) 

PRZBF(3) 

PZETBF(3) 

PRZBP(3) 

PFBLP 


[tU9] [TYB] [T51]|o| 


The aercxlynamic force vector is 



Mode shape data is available in array fora* 



where the array dimension definition is 
(2, 3, 20) (COORD, MODE, STATION) 

MCOEF (2, 3, 20) 

MCOEFP ^2, 3, 20) Prime values (derivatives) 

FBIB (2, 3) Feather bearing inboard 

FE03 (2, 3) Feather bearing outboard 

5. 2. 1.1 Rotor Blade Spatial Integration . - All force and mass integrals over 
the blade are computed from the inboard to the outboard stations. The numer- 
ical algorithm used is trapezoidal integration. Let FI stand for the function 
(integrand) to be integrated emd F for the integral. 

Then the algorithm proceeds as follows. 

• Evaluate FI^ i.e. FI(X^) 

• Set F^ i.e. F(X^) = C 

• Integrate from k=2, , 

/ k _ k-l\ 

F^ = F^“^ ) [fi^ + FI^"^] 

This process is pict\ired in Figure 22. 

The integral equation presented above for force and mass terms are 
straight. However, the current subprogram treats each function separately. 
Therefore the correspondence between program symbol and function is 
presented. There €ire currently 62 functions defined for each blade. 

The array size is F(66, 7). the integrands for each function are 
defined in the FI array. The definitions of F for each blade are 
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aerodynamic loading only 


Mass Integrals follow. 


oA ) aerodynamic loading only 
BL I 


^6 = 


BLE n 


^19 



^20 

= 

”Vi 

F 

21 


”V2 

F 

22 

= 


F 

23 

S 

\k 

*11 

C\J 

s 

”*2*1 

F 

25 

= 

“*2*2 
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Continuing Mass Integrals: 


^26 


M . 

^27 

= 

”a A 
*3 2 

^28 

s 

\ A 
3 3 

^29 

= 

M, 

’*'bla^ 

^30 

s 

BLAg 

^31 


%LA 

^32 

s 

«» 

BL f 

^33 

s 



* 

\u. 

^35 

s 

%U3 

^36 

= 

"^BLA, 

^37 

= 

BLA^ 

^38 

= 

"^BU, 

^39 

= 

BL0^ 

^0 

~ 

”z 

^BLA^ 


s 

“z 

bla 

^2 

= 

”z 

BLA 

"U3 

s 

^Z 

BL^.^ 

"uu 

= 

’bl 1 

^5 

= 

A 

%L 2 

^6 

= 

A 

*BL^3 
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^7 


% ♦ 
^BL^f 

^8 

a 

®BIA^ 


a 

®ELA 

"50 

a 

Mg 

^BLA^ 

"51 

= 

BL’^f 

"52 

= 

M. , 

li th 

BL^BL 

"53 

= 

®bl'^bl 

"51* 

s 

My ^ 

BL BL 

"55 

s 

^BL*^BL 

"56 

= 

^BL®BL 

"57 

a 

M* . 

Vl 

"58 

= 

M* , 

V2 

"59 

= 

M* . 

V3 

"60 

a 

M, , 

T T 

"61 

= 


"62 

a 

V*f 


5.2.2 XTEFMO . - This is a specialized matrix routine, subordinate to SWEEP. 
Within SWEEP, the equations for 

• Mbl • Mb. 

n n ' n 

require the assembly of complex matrix terms which are made up of transforms 
of three angles z' , and Subprogram XTERMO and its multiple entries 

provide the calculation of these terms. The program code has been computer 
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generated by an application of IBM's FORMAC program where the number of cal- 
culations is luinimized. The subprogram entries are described below. 

XTERMO - an initialization entry where trignometric terms of the three 
transformation angles are computed. 

The three basic transformation are, in argument order; (Zp^, ^'p^, M 



Z'^ 

0 

sin Z'^ 

[t ' ] = 

^ ZFA^ 

0 

1 

0 


-sin Z'^ 

0 

cos Z'^ 


cos Y'^ 

sin 

Y' 0 

FA 

[t ' 1 = 

^ YFA^ 

-sin Y'^ 

cos 

Y' 0 

FA 


0 


0 1 




1 

0 

0 


:os 


-sin 


sin 

cos A(f>^ 


The transforms computed on the various entries are 


XTERMl 


[T] = ^'^ZFA^ 


XTERM2 


Ft! = Ft* Ft' Ft Ft' 1 Ft* i 

Lij L^zfa^ ^^YFA^ ^ ^ YFA^ ^^ZFA' 


The dot over the matrix indicates differentiation with respect to the angle. 
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[t1 = Tt' Tt' fT Tt* 1 (t* 1 

LiJ ^ A(J)f^ '■ YFA^ ^ ZrA-* 


XTERMi* 


XTERM5 


XTERM6 


XTERM7 


[T] = Ft' fT' fT Ft* 1 Ft' 1 

L-J Li2FAJ ^^YFA^ ^ A(}ifJ ^^YFA^ ^ ZFA^ 


^"^ZFA^ ^"^YFA^ ^"^A<{)f^ ^^ZFA^ 


Ft] = Ft' Fi" ]” Ft r Ft' 1 Ft' 1 

LiJ ^ YFA^ ^ A4if^ ^ YFA^ ^^ZFA^ 


■*■ ^'^ZFA^ ^"^YFA^ ^"^YFA^ ^"^ZFA^ 


FtI = Ft' l'^ Fi" l'^ 'i Ft' 1 Ft' 1 

^ ^ ^ ZFA^ ^ YFA^ ‘ AV^ ^ YFA^ ^-^ZFA^ 


+ FT' ]'^ Ft' l'^ Fi Fi* i Ft' 1 

^ ZFA^ ^^YFA^ ^ A4,f^ ^^YFA^ ^ ZFA^ 


Ft] = Ft' Ft' Ft i Ft' 1 Ft' 1 

LiJ l^zFA^ ^^YFA^ ^ A1)^^ ^^YFA-' ^^ZFA^ 


+ 2. fi' Ft* f'^ F"" Ft’ 1 Fi* 1 
^ ZFA^ ^ YFA^ A<{.f^ ^ YFA^ “-^ZFA^ 

+ [Ti„ ]'^ Ft* ]'^ Ft l'^ Ft” ] Ft* 1 

ZF>. ^ YFA-' ^ ^‘YFA^ '■ ZFA' 
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XTtSMS 


XTEIg{9 


XTZRMA 


XTERMC 


IT] = [T‘p^] 

'*' ^"^ZFA^ ^"^ZFA^ 


[T] = [T^p^] [T^^] ^"^YFA^ ^*^ZFA^ 


'*' ^ ^"^ZFA^ ^"^ZFA^ 


* ^"^ZFA^ ^'^YFA^ ^"^ZPA^ 


[T] * [T'p^f [T-p^] 

■*■ ^*ZPA^ ^*YFA^ ^^YFA^ ^*^ZFA^ 

^"^ZFA^ ^^'b’A^ ^"^YFA^ ^"^ZFA^ 

^"ZFA^^ ^"^YFA^ ^‘YFA^ ^^ZFA^ 


t%A^ 


XTERHD 



XTERME 


[T] 


[T* 1^ [T* 1^ 

^ ZFA^ ^ YFA^ 


5.2.3 COULOM. - Coulcnb friction is calculated toy sutoprogran COULOM. COULOM 
is used for svashplate and feather bearing friction. 

The output function is of the form shown in Figtire 23. 



Figure 23. - Coulomb friction. 


The routine is executed by 


CALL COULOM (F, FC, RL, RT) 


F - output friction 
FC - friction limit 
RT - argument such of 
RL - argument limit 

5.2.U Main rotor blade aerodynamics . - The following three routines, CLCDCM, 
STALL, and XTRPl provide main rotor aerodynamic coefficient calculations. 



5.2.14.1 CLCDCM. - This routine calculates the aerodynanic coefficients 

=D * «> 

by linear interpolation. Up to two sets of data tables are available in 
labeled COMMON/TAB/. This routine is derived from a similar routine in the 
Army's C-8l computer program. 

Arguments List 

IDT Data table set number. 1 or 2 

ALPHLD Angle of attack for CL and CD curves. Urdts are degrees 

ALPHCM Angle of attack for CM tables 

MACH Mach number 

CL 
CD 
CM 

5.2.U.2 STALL . - Subprogram STALL applies the dynamic stall model as pre- 
sented in section 6. 2. 3.1+ of Volume I. These equations are distilled in 
Figure 2U. 

Arguments List 

ISET Aero table set nvunber 

A Angle of attack, rad. 

M Mach number 

AYAW A radial flow angle. Tan ^ (Ug/U^) 

K Blade element reduced frequency 

K - Ci/2V 

10 
















Arguments List 


CLR 

CDR 

CMR 


Output aero coefficients 


5. 2.1*. 3 XTRPI . - This is a special linear interpolation routine designed for 
speed. The interpolation technique deos not require the classic table search 
to find the data of interest. Ihis routine does require that the tabulated 
function be evaluated at constant argument intervals. 


Algorithm - If x is the first argument value and the data is evaluated 
at Ax intervals, then the location within the table of bracketing func- 
tion values for a value of the argument x is n where 


n = 



The brackets indicate the largest integer whose magnitude does not exceed 
the magnitude of (x - x^)/ax. 

The interpolating formula, based on the point-slope formula, is given 
as 


y = y + dx 
n 



y 


) 


where 


dx = 



- n 


This is shown to be correct if 


x-x *n'Ax+{x- x„) 
o “ 

is substituted in the equation for dx. The quantities x and Lx are stored 
with the table of fiuiction values. The table, T, must be of the form 
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T(l) 

X 

o 

T(2) 

Number of x breakpoints - 1 

T(3) 

1/Ax 

T(U), ... 

Table data of increasing x. 


5.3 BLADSF 

Subprogram BLADSF computes the geueralized force vector for the blade sub- 
system. The single argument to the routine, DELQBS, is that force vector. 

The blade subsystem force vector is defined 



where for each blade n 



for all modes m = 1, .... 


Consider first, the bending mode equation. The vector has been 

computed in SWEEP where ^ ' 
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(eleaents of the array F{66. 7)1 



The potential energy term is given by 



Tvo damping terms are computed. The structural damping is given by 



The proportionality factor , C^, is controlled via input. The spring constants, 
, are also input. 

A lag damper with relief valve is also modeled. For each blade n 



Now if 




then 


^LAGl ^n 

c 
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otherwise 


«LAG = "lag, ^3 " "lag, (^n ‘ ^b) 

1 2 ' c ' 


where 








Finally, the damping term is 


"CD (S') 


In the above equations, the following parameters are input 


3A 


m 


m = 1, 2, 3 


YP(1) 

YP(2) 

YP(3) 


= YDPI 

"lagi = 

"lAG2 = "^"2 

The fourth mode in the set cau be pitch horn bending or dynamic torsion. If 
Ngp y 0 then 
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If ^0 then 
T 


AQ =F — K •li — C. ‘i 

“^3^ 17. n V ®T "'i- 

n ’ T n T n 


5.1* HUBF 


Subprogram HUBF computes the force vector relative to the hub, DELQH. 

The hub is not a generalized coordinate ?ab£>scem, but does provide e transi- 
tion of forces from the main rotor to other subsystems. Equations of inter- 
est follow. 




0 -r 
r 

-q 

0 -r 
r 0 

-q P 




— 

/ V 

q 


3x3 from 

(" 

-p 


upper right 
quad of 


0 



M 

— 

H 

L 

\ / 

q~ 


3x3 from 

(") 

-p 


lower right 
quad of 

H 

0 



n 


H 

— — 

\ 1 


lAF| 


SP., 
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and 
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The main term in the hub force equation is the main rotor force relative to 
the hub. 



The blade generalized forces are computed in SWEEP as elements of the F array. 
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5.5 SWASHF 


Subprogram SWASHF calculates the generalized force vector for the swash- 
plate subsystem, DELQSP. The swashplate subsystem consists of a maximum of 
three coordinates. 


qgp = (4' t 9 . Z)gp 


The force equation is given as 



is the feathering moment computed in SWEEP. 


The frictional force vector 



''f(^sp^ 


:f(6sp) 


where the frictional functions are determined by entries to the subprogram 
COULOH. 

The term |af| as written above is a (3 x l) vector with the stated 
coordinate ordering. The definition is given below in the standard (6 x l) 
form for clarity. 


Ill* 




AF 


SP 


0 0 0 1 0 0 

0 0 0 0 1 0 

Lo 0 1 0 0 0, 



(cose (I^p + qr (I22 - lyy)) + sine • * r)gp^ 


= {(%<! + • (Ij^- l22))sp 


(m Z) 


SP 


The potential energy and dissipation functions are presented below. The 
angular coordinates are treated separately from Zj,_. 

O i 
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Let 


[Kgpl 


2 2 
K. CCS ii + K„ ;jin >P 

*SP "SP 


K . - K„ sin ^ cos i‘ 

*SP ®SP 


/K. - K. \ Sin^ 

I ♦sp ®sp) 


Kw sin^ 

^SP SP 


and 


l=spl = 


'■ K ■ <=6 ) 

V ^SP SP/ 




cos 4i sin \p 
c c 


-sin cos 'J' 

c c J 


Define a control load 



-%s ■ 


K, • V 


cs 


cs 


Nov 


| au 

|3q 


SP, 




(*<’sp) 

0 


SP 


■ K] 


(h ( 

|3q J 


tCgp] 



A stop load is modeled. 


Il6 


CCS >P 

c 

***’”**’ 

cos“ 

c 



-<■ (»SP^ * * 's, s? 


= L • 


|*SP| 
I ®SP 1 


^•■'s.sp [(y-y)"''- /(y* y) 


otherwise 


r=i 


I su I _ tau I 

I ) SP I I SPj^ 


IJ 

risp 


- ~ 7 * "CP * 


otherwise 


‘'SP SP 


^SP ■ ^ISP ) * 
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where 


^ISP “ * ^G1 


_3B 

32 


= C 


SP 


^p 


SP 


- K 


['...] 


SP 


*SP ^SP 


- R„ 


['2.2] 


SP 


"SP 


SP 


The above eqiiations contain a great deal of data. The equation symbols, corre- 
sponding program symbols and input addresses are listed. 


Symbol 

FORTRAIl 

RA 

^c 

CHI 

119 

^CS 

QKXCS 

123 

*SfCS 

QKYCS 

I 2 U 

ir . 

^SP 

KPHCOH 

376 

^^SP 

KTHCON 

377 

c. 

CPHCON 

378 

■"SP 

Ca 

CTHCON 

379 

“SP 

^s,SP 

GASTOP 

1276 

^s,SF 

GKSTOP 

1277 

*SP,BR ) 

RLG 

116 

®SP,BR ; 

”fr* ^sp 

( FCG 


117 

”fr’ ®sp 

) 


K 

QKGZl 

137 

IZSP 

z. __ 

ZGl 

lUl 

ISP 

^ 2 Z 

QKGZ2 

iko 

F 

FIDDLE 

1 U 9 U 

c 


QCGZ 

138 

“SP 

R 7 » R 7 

DGDHG 

1263 
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The reader can reference Section 5.9 of Volume I of this report for addi- 
ticnal control system explanation. 


5.6 SHAFT!' 

Subprogram SHAFTF computes the generalized force vector, DELQS. The force 
eq-iation is 


fe1 hk N -N Ms 

The hub force vector is computed in HUEF and passed as an input argucient. Since 
the shaft coordinate set can be a subset of the standard six coordinates, the 
Kg and Cg coefficient matrices have been asseabled in the data preparation 
subprogram DTPRSP. The matrices are Ng x Ng symmetric. The reader can refer- 
ence Section 5.7 of Voliane I for further information on the shaft force 
equations . 


5.7 REFF 

Subprogram REFF computes the reference system generalized force vector, 
DELQF. The force equation is 







The hub force vector is computed in HUBF and passed to RLFF via the 
argument list. The fuselage force vector is given by 



* 0 -r q* 
r 0 -p 


(3 X 3) 
upper right 
quad of 


V 

q 

_-q p 0_ 

f 

L«of J 


_r_ 


’ 0 -r q 


(3 X 3) 

j 

1 ") 



lower right 



r 0 -p 


qixad of 





M „ 

: 


.-q P 0 . 

F 

.OF J 


^ r ) 


F 


The fuseleige aero loads depicted by are computed in subprogram LOADS. 

The tail rotor cmd engine contributions are 




’’ref 
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where 


r 


R 


= V 


R 


r 


R 


R 


The above equations can be referenced in Section 5-8 of Volume I. 


5.8 LOADS 


Subprogram LOADS computes a total vector of nonmain rotor aerodynami 
loads in fuselage coordinates. This force vector is the primary vector o 
the reference system force vector. The loads summation is given by 


W,) 



where { 1 * 3 }^ stands for total body loads, {1^} is tail rotor loads, and -.Lp}^ 
is the thrustor T ads vector. Program symbols are 



QL0ADS(6) 

FN(6) 

FTR(6) 

FP(6) 


SI(^ (density ratio) 


The equations are straight forward but extensive. The reader should refer to 
Sections 6.3, 6 .U, 6.5, and 6.6 of Volume I for equation development. 
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5.9 ROTORF 


Subprogran ROTORF computes the generalized force for the rotor (engine) 
subsystem, DELQR. The system is a one-coordinate system and the equation is 
given by 





3'*’ ENG 
R 








where 








My = 0 

^NG 



My = My 

^ Nr 


% 


o 


M 


$ 


f 



M. 


f 

n 


The reader can refer to Section 5.10 of Volume I for equation development. 
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6. MASSSP - LEVEL 1* EXECUTIVE 


MASSSP is a level U executive control program. Its prinary output is 
the generalized mass matrix. The elements of MASSSP are depicted in Fig- 
ure 25 . As in the development of the force vector, a subsystem approach is 
followed. 

The generalized mass matrix is symmetric positive definite. Thus, only 
half the matrix is computed. A schematic of a complete mass matrix is pre- 
sented in Figure 26. 

This subprogram controls the assembly of the generalized mass matrix. 
Each BKsdel subsystem is developed individually as required. The required 
subsystem coupling matrices are stored in the generalized mass matrix via a 
utility routine, RMSTOR. The subsystem diagonal matrices are stored by the 
utility routine, DffiTOR. The blade subsystem mass matrix is handled directly 
by its development routine, BLADSM. 

The subsystem matrices, their sizes, and the location of each within the 
generalized mass matrix are listed below. 


Subsystem 

Matrix Size 
n, m 

Location Within Gen. 
Matrix, M(i, j) 

i 

tj 

Swashplate 

Mspbl N3S) 

LSP 

1 

Shaft 

»SBL 

LS 

T 

a. 


«SSP 

LS 

LSP 

Reference 

’'HEFBL 

LREF 

1 


•Wsp 

LREF 

LSP 


(miT, 8S) 

LREF 

LS 

Rotor 

''m f™. 

LR 

1 


Mp2p (NR, NSP) 

LR 

LSP 


Mpg (NR, NS) 

LR 

LS 


Mppgp (NR, NREF) 

LR 

LREF 







f'O 

fr 



MATRIX ROUTINES 


MXMULF 


MTMULF 


CONMfS 


Figure 25. MASSSP - level U hierarchleval chart 















Fig^ire 26. 




The flow of MASSSF is depicted in Fig\ire 27. Routines for each model 
subsystem are easily recognized. Two other routines identified are ZERORM 
and HUBM. These routines compute intermediate mass matrices which are used in 
developing the required coupling matrices. 


6.1 DMSTOR 

This routine stores a subsystem diagonal mass matrix into the generalized 
mass matrix. The subsystem and genersilized mass matrix are in symmetric stor- 
age mode form. 

Arguments List 

L Subsystem coordinate vector location within the generalized 

coord, vector 

N Subsystan size 

MSUB Subsystem diagonal mass matrix 

MSP Generalized mass matrix. 

Subroutine DMSTOR is depicted in Figure 28. 


6.2 RMSTOR 


This subprogram stores a subsystem rectang\ilar coupling matrix 
into the proper locations within the generalized mass matrix which is^in 
symmetric storage mode. 

Argxments List 


II 

M 

UMAX 

MSUB 

I 

J 

MSP 


Submatrix dimensions M 3 ijg(N, M) 

Row definition in DIMENSION statement for MSUB 
Submatrix 

/ Logical location of MSUB(l, 1) within the generalized mass 
' matrix 

Generalized mass matrix. 
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*^BL 




**SBL~*”**^ 


Mssp— 




M, 


REFBL 


— M®* 


M, 


REFSP 


^M®* 


^REFS 


M®* 


12T 









Figure 27. - Concluded. 
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Figure 2«. - DMSTOR flow diagram. 
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The problem ^s, given the size of the rectangular submatrix {n, m), to 
locate the origin within the logical generalized matrix, Then 

Mgyg (k, h) - M (i - 1 ♦ k. J - 1 + k) 

Ic “ •••) n 

~ Xy •••y ni» 

The real mass matrix is in symmetric storage mode, i.e., in vector form. How- 
ever, an element, S, of the vector corresponds to the matrix element qr where 

s = (q) (q - l)/2 + r 

The final algorithm is then 

Mgyj (1^. h) - V(t - 1 + h) 

where 

t=(i-l-tk-l) (i-l + k)/2 + J 

for 

k*~l, •••, n 

h *— 1, m 

This is depicted in the Figure 29. 

6.3 BLADSM 

The subprogram BLADSM assembles the elements of the blade subsystem gen- 
eralized mass. The use of intermediate matrices is minimized because of the 
specialized form of the subsystem. Data is stored directly into tne general- 
ized mass matrix, MSP, which is the only quantity in the argument list. 
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The mass matrix for the total blade system can be depicted as 



The symetric mass matrix for each blade, n, can be further expanded as 



The subscript. A, identifies blade bending modes. The number of bending 
modes, is optional. The addition of another mode, S» is also optional. 

If the input NBP = 1, then £= where Bpj^ stands for pitch horn bending 

(elastic feathering). If the input NPT = 1, "then 6= where stands for 

dynamic torsion. Also, note, that the number of blades, K^, is optional. 

Current limits on blade coordinates are 

NB T blades 
b 

N NMB 3 bending modes per blade 
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«BP 


••4,7 


KBP 1 elastic feathering mode per blade 
NPT 1 dynamic torsion mode per blade. 


Ngp and N are mutually exclusive. 


!?T 


The mass matrix definitions to be given below are in terms of mass inte- 

greuls developed in subprogram SWEEP. Blade force and mass integrals are stored 

in an array F. . Expanded definitions can be found in the section on SWEEP. 

1 ^ n 

For each blade, n, the elements of the symmetric matrix [Kg. ] in sj-rmetric 
storage mode are 


% ' ^22i-J,n 


J = i. 


‘ mb 


and if N^p 0 



or if S* 0 

" ^56+J,n 
^33 ~ ^60, n 

otherwise 



0 


Mgg z 0 


Qp 

The utility routine, D^B'POR, is used to load each [Mg^ ] into [M]‘" as each 
blade is evaluated. 


133 



6.1* ZERORM 


The subprogram ZERORM computes a set of intermediate matrices for use 
other routines. These include and Mq™.. The rationale is to collect 
blade-dependent data in a pseudo syst«n wni^ is called zero of the rotor, 
the subscript (OR). The matrix definitions are 


w-s M. fe! 

n=l ^ ^ 

(6x6) 


where 





1 



- 



0 

0 1 
1 

0 

^6 

■^55 


0 

”bl 

0 

■^56 

0 

^ 51 * 


0 

0 

«BLi 

''ss 


0 

0 

~"76~ 

1 

'’’”1 

^XXBL 

0 



^56 

0 

-" 51*1 

0 

^YYBL 

P 

■53 


-"55 


1 

0 1 
1 

^52 


^ZZBL 


and 


where 



by 

all 

thus 
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Continuing with the definitions. 


BLAm I “oBLg J ^ 


^ere for all blade beading modes, m. 



If # 0 then 
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Otherwise 



is not defined. 


The partial derivatire 



is of the special fora 



Therefore, 

Re-define 


the iiQ>leaentation of and ] proceeds 


as follows. 


"b 

[”or]" S [**OR]n 


where 



r I n 

^11 [ ^12 

^ 21 ' ^22 

- I J 


and for each blade, n. 


■"iB ' [Vbl]^ ["bLOjJ ['^R-Bl]„ 
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22 


' 





0 








n 


Finally, partition 


[”0BLBL j ^ 


C, (3 X N ) 
X m 


C_ (3 X IJ ) 

d m 


then 



Note that the subscripted F's are mass integrals developed in subprogram 


SWEEP. The blade mass, inertias I 




"Ti 


BL 


idLt. 


are calculated 


BL 
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in subprogram INIT as well as the transformation matrices from rotor to blade, 
[Tp and the partial derivative 


6.5 HUBM 


The main rotor hub serves as another pseudo system for gathering informa- 


tion used by other subsystems. Subprogram HUBM accepts two matrices as input 
directly through the calling sequence [Mor] and fMnPRT.l. It computes and 


outputs three matrices ^**ORH^* definitions of these 

matrices follow. 


[«orbl] 

[“a] ["ora]* [’^]h* 

["ora] ' ["o"] {5^ I 

The transmission mount mass matrix, is constant and computed once 

in subprogram INIT. The mass of the swashplate in hub coordinates is 
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The partial /3Tg| coupling the rotor 

prograa PDERIV. Other definitions include 


to the hub is computed in sub- 


OSP 


7 ' 

OSP 


m. 


SP 


^P 


YY, 


SP 


ZZ 


SP 


^OSP * ^SP 
ZZSP 5 RA(11i69) 

SPMASS = RA(139) 

IXXSP = ra(36i) 

lYYSP = IXXSP 

IZZSP = RA(118) 


6.6 SWASHM 


This routine con^jutes the necessary swashplate mass matrices. These in- 
clude the svashplate-blade system coupling matrix and the swashplate-swashplate 
symmetric mass matrix. The coupling matrices fWopgp] 
puted. The matrix definitions follow. 



1 ”SPBL 

I Nd 
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where 



If Ngp ?^0 (pitch-horn bending) 


M 

X 




4> 


f 


or if N. #0 (dynamic torsion) 
T 


M = M. * 

X y $ 

f T 


otherwise 


M = 0 

X 

Further definitions include for each bleide n 


” Vm " ^19+m,n 


»®~1* 


mb 


M, 


= F 


f f 


19, n 


lUO 


^'^’f'J>T " ^6l,n 



The following partial derivatives are also required 






the partial of feathering angle wrt the generalized 
swashplate coordinates 


If Ngp = 1. 


If Ngp = 2. 


If Ngp 3, 



3<t^f 

30SP 


39sp 


> 


94*^ 

3Zgp 


Alsc, the partial of the feathering angle wrt the generalized rotor coordinate 
is 



0 . 



where 


3^j. 

= A^g sin '*'f) ■ ®1S ^r) 
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The swashplate - swashplate matrix is 


^WP ~ * ^^P Q 




The swashplate coordinate set is not standard nor are they grouped in a stand- 
ard order. With this in mind, the following mass matrix and partial derivative 
matrix are based on the coordinate definitions: 




lk2 



1 


The program implementation of the above definitions follows. If "SP = 
i.e., Zgp coordinate only* then 


N] 0 “ 

(1 X 1) 


[^PBs] S 

(1 X 1) n=l 




n 


otherwise in .‘Symmetric storage mode 





(1) 

(1) 

(2) 

(3) 

(2) 

(3) 

> 

(U) 

(5) (6) 

ih) 



(5) 



(6) 


^XX„ ®SP ^ZZ-p ^SP 


'^^SP - Xp 


0 

0 


m 


SP 


In symmetric storage mode, Mspgs ® vector of length Ngp ‘ (Ug^ + l)/2 

W^>-I M 

J ” Xy • * • % ^SP 

lt"“Xy eeay J 

L = J • (J - l)/2 + k 





Finally 


[«swpj__= [«SPo] * [”sPBs]^ 

m 

0*1 “sp^^sp 

The remaining output matrices are 


' ksp ] * [““"=''1 



The partieds are not defined explicitly. The particular form of [M^5pl 
will be determined as follows. 



If NSP = 1 then 


= ( I j t I ^ 


["hSp] ' j j 


I '’S'’ 1 


0 

0 

“sp 

0 

0 

0 


otherwise define 


[*^sp] 

(6 X 3) 


0 

0 

0 


^XXsp ''hsp^^ ®SP " ^ZZgp ■'hsp3^ 


0 

0 

0 


sin 3gp ^YY '^HSP21 




0 

0 

SP 

0 


cos Sot, ♦ Ir 


■•XXgp HSP,2 sp HSP 


SP ‘-32 ““'sP ^YYgp^HSP22 


^XXgp "^HSF^^ ®SP * ^ZZgp *^HSP33 ^SP ^YY^p "^HSPSS ° 


and the final 


[^sp] “ [”hsp] 


6 X N, 


SP 


6 X N, 


SP 


The matrix [MoBL^f^ defined in ZERORM. The partial |3 "o?L^^'OR( is 
not available explicitly but is defined as 
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6.7 SHAFTM 

CuV;rogrr- SHAFTM conputes the mass matrices required to represent the 
shaft as a sciJel subsystem. Matrices input include 
[My], Matrices output are defined below. 



r 1 

[«SBl] - j 5Tj j 

[•Wj 

[«ssp] = j a.f 1 

>T 

1 ["“S''] 

K T ( ) 

T 

r 1 


["-] 

[>S,s] = [«h]' 

1 ’’OH i 
j-s ) 

[”0Bs] ' 

1 


The partial defined in PDERIV. 
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6.8 REFX 


Subprogran REFM cc»putes the mass matrices necessary to describe the 
reference subsystem. Input matrices include [%bl1» 

[Mg], and [)^] . Output matrices are defined belov. 



^*Ws j 

[•w] 





I 

I 


1U7 



The matrix [Mo]p is a constant and ccaiputed once in INIT. I^ie partial 
jdTOH /^Tp£p| is computed in PDERIV. The partial connecting the fuselage to 
the reference systen is defined as 



6.9 ROTORM 

Subprogram ROTORM ccxaputes the mass matrices necessary to describe the 
rotor subsystem. The subsystem matrices are defined below followed by im- 
plementation information. 


where 



1U8 





■ 1 ’’s 1 1 "°“®] * 

m 

( ) 
i 


•|"”t w 

1 3^0SP \ 

\ ^^REF ) 


W' 

1 ^^‘ENG ) 
( ^^REF ) 

■l?l hlfc! 

1 4" ( ['.] 1 4" 


■| [ 1 I 


+ 




♦ 



The following matrices are input, tMopspl » iMorsK t^RREFl» 

[Mor]* Other matrices and partial derivatives are defined below. Finally, 
the implementation of the individual terms in the above equations are pre- 
sented. Partials of interest include 




1 49 
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Terms are now defined ais 



Note the forth element is approximately Zero. 
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cos $ p sin egp cos 9gp 


3t 

ENG 


["bep] ■ {■ 


, 0. 0 I 


^ENG ^XX 


. 0 . 0 ^ 


ENG 


HISl' I"' 


I iC"! ["“=] i 


ENG XXjjjg 


TR 


I22 ^SP ®SP 



0 , 0 , 0 , 0 . 0 , 


- ^SP ^SP( 


^"^OSP I i 


at 


3T 


RF 


0, 0. 0, 0, 0, - 


^ZZ ^SP '^SP ®os ®os ®s 

SP 
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7. MATRIX ROUTINES 


Throughout the documentation, reference is made to Full Storage Mode and 
Symmetric Storage Mode. This refers to the manner in \rtiich a matrix resides 
in core storage. A matrix stored in Full Storage Mode assumes the normal 
storeige attributes as dictated by FORTRAN for two-dimensional arrays. 


Symmetric Storage Mode is used in order to minimize memory requirements. 
Only the elements on and below the main diagonal of symmetric matrices are 
stored. The order and occurence of these elements in core are as follows: 
Assume [a] is a symmetric matrix and [B] is a vector. Then 


[A] 


““ — 


^ — 

^11 ®12 ®13 


hi 



♦ 

®21 ®22 ®23 

— 

®21 ^22 

^31 ^32 ^33 


^31 ®32 *33 

— 


— ^ 



22 


®31 

*^32 


[B] 


An n by n matrix in Symmetric Storage Mode is reduced to a vector of length 
n(n + l)/2 where the element ij can be found as the element k of the vector B, 
where k = (i (i - l)/2) + J for i - J. Storing a matrix in this manner effects 
a savings of n (n - l)/2 core locations. 

The Cholesky decompositicr routine, CHOSKY, is specifically designed to 
operate on a matrix in Symmetric Storage Mode as is the output routine, 

MXOUT. 


The routine, CO?iMFS, will convert a matrix from Full to Symmetric Storage 
Mode. The remaining routines, MXMULF and MTMULF are standard matrix multiply 
routines. 


7.1 CHOSKY 


Within REXOR II, the equations of motion are stated 



where 


|Aq[ = [M]"^ |AF| 


M and AF are given. Aq can be ccanputed by first inverting M then performing 
the indicated multiplication. A more efficient method is to solve for the 
product directly by solving the linear system 


[M] {Aq} = |AF| 


for |a^|* It is further known that the mass matrix, [M], is positive definite 
and symmetric. Subprogram CHOSKY is a general algorithm for the solution of 
simultaneous equations of the form 


[A] |X| = |b| 


where the coefficient matrix is positive definite, symmetric. The algorithm 
is a Cholesky decomposition of [A]; followed by a forward-backward substitution. 
The algorithm is presented below. 

Cholesky method for symmetric, positive definite matrices . 

Theorem: Let A be symmetric, positive definite. Then A can be factored in 

the form 


T 

LL = A 


where L is a lower triangular matrix {i.e., L = (l^j) where 1^^ = 0 for J > i), 
Cholesky Method; Let A be n x n, symmetric, positive definite 


A = (ajj) 




Assume A is factorable 


A = LU 
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Then 


“ij ' 2 *Ik \j 


“ - "W ' ‘jk 

then 


a 


ij 



^Jk 


I 1. ** Xy •••^ n^ J ** X^ n 


r *•- L is lower triangular which impXies = 0 for k > i. Therefore» 


a 


iJ 



*'jk 


X * Xy n^^ J “ Xy ■••y n 


This equation forms *he basis of the decomposition. The eXements of L are 
found as foXXows. 



aXso 


®XJ ' ^XX ^JX 


leads to 


^JX ^XJ^^Xl 


»J“2, •••1 n 
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and 


a 


ii 


•z 


i-1 



k=l 


• 1~2| n 


can be solved for c . . . Therefore 

ii 



k=l 


Finally, 




i-1 

=E 

k=l 


“•ik *-Jk ^ii *-Ji 


, i * 2, n 


giving 


£ 


Ji 



,j=i+l, •••, n 


With the decomposition of A accomplished, the system of interest can be 
solved as follows. 

Given 


A X = b 


Substitute 


A = 
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giving 


T 

LL X = b 


Define 


T 

g = L X 


Solve Lg = b by forward substitution. 
Namely 


®1 ^I'^^’ll 


and 


«i 



i ~ ...f n 


Finally, X can be found by backward substitution. 



The above equations form the algorithm. 
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The arguaents in the CHCSKY subprogram calling sequence aie 


Argmaents List 


A An Input vector of length l)/2 containing the i; x N 

I>ositiv€ definite symetric matrix stored in symetric storage 
mode form. 

Input matrix of dimension K x L containing the right hand 
sides of the equation AX = B. 

On output, the K x I solution matrix (X) replaces (B). 

Order of [a] and the number of rows of B. 

Row definition of B in calling program. 

Humber of right hand sides (columns in 3). 

Error indicator 

= 0 normal return 

> 0 indicates the input isatrix is algorithmically 
not positive definite. lER with contain the 
index of the row which failed. 

The positive definite test is as follows. While forming the eleaents of 
the decnuposed matrix, L, the argument of the square root f-jnction is tested, 
i-f 


B 

H 

M 

L 

lEP 


i-1 



then the procedure is terminated. The coetficient matrix. A, is destroyed 
on exit. 

The subprogram CKOSKY is c«nputed in double precision on IBM hardware. 
This means the matrices A and B are double precision. The CDC version is 
single precision. 



7.2 CONMFS 


This routine is a matrix storage mode conversion routine. It converts a 
full symmetric matrix to a matrix in symmetric storage mode. 

Arguments List 

A Input matrix of dimension N x K. 

N Order of matrix A. 

lA Row dimension of A as specified in the calling program. 

B Output vector of dimension H • (H ♦ l)/2 containing 

matrix A in symnetric storage mode. 

An n X n matrix 5.n symmetric storage mode is reduced to a vector of 
length n(n l)/2 where the eleaent ij can be found as the element k of the 
vector B, where k = (i (i - l)/2) + J for i >J. 


7.3 -HXOUT 

This routine prints a matrix which is stored in symmetric storage mode. 
Arguments List 

A Name of matrix to be printed. 

M Number of rows and columns of A. 

K Number of lines already printed. 

Output is to FORTRAN unit 6. The routine will page eject after 55 lines. 

7.U MXMULF 

MXMJLF is a matrix multiplication routine, vbere the matrices arc in F\iU 
Storage )tode. 


[c] = [r 


Arguments List 
A 
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B 


(LXM) matrix 
(MXN) matrix 



L 


Maximm value of first subscript of A 


M Maximum value of second subscript of A and first subscript 

of B 

N Maximum value of second subscript of B 

lA Value of first index in DIMENSION statement for A 

IB Value of first index in DIMENSION statement for B 

C Resultant matrix, (IXN) 

IC Value of first index in DIMENSION statmaent for C 

The multiplication is performed in double precision. 


T.5 MIWULF 

This subprogram performs a matrix multiplication of the transpose of a 
matrix A by a matrix E. The matrices are in Full Storage Mode. 


[C] = [A]*^ [B] 


Argum«its List 

A (LXN) input matrix 

B (LXN) input matrix 

L Number of rows in A and B 

M Number of columns in A and rows in C 

N Number of columns in B and C 

lA Row dimension of A as specified in calling program 

16 Row dimension of B 

C Resultant matrix (MXN) 

IC Row dimension of C in calling program 

ihe multiplication is perforaed in double precision. 
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3. program hODIFICATIOH HIBTS 


“Ihis section is not intended to be a coa^jrehensive progranaaing guide, 
but rather a few hints concerning the extension of the prograa in the area 
of input /output. 


8.1 Adding Paraoeters to Input Set 

All prograa data which is user changeable via input resides in a data 
array, RA, which is assigned to the labeled COMMOS/IHFUT/ as follows 

CC»9«OS/INPUT/RA(3000), TITLES (15,1*) 

Not all of the 3000 addresses have been assigned to input. 'Hiis can be seen 
in Table 3-1 of Volwse III, where all addresses are defined. Those which are 
undefined as of this writing are indicated as "open". If one wishes to 
define a new input, -^ind an undefined address in the table and equivalence the 
new parameter to the appropriate address. For example, assume a new input 
is required, X5EM, and am open address is found, 20. Hioq sinply declare 

EQUIVALENCE (RA{20), XHEW) 

in those routines necessaury. A copy of the labeled C0l®*0S statoaent 

COMMON/ INPUT/RA( 3000), TITLES (15,1*) 

is also required. Another fact should be pointed out. The aurray, RA., is a 
REAL variable and all inputs are read as REAL data. If am input is required 
in integer form, then define the new parameter across an equal sign, 

IP = RA{20) 

Accurate records of assigned addresses should be kept to avoid conflicts. 


8.2 Adding Variables to Trim Set 

Program variables available for trim control or function definition are 
found in the table which drives the RETRVE/STORE routine. The table is a 
form of indirect addressing. The table can be easily extended. Currently, 
variables defined in blank COMMON and /INPUT/ labeled COMMON can be accom- 
modated by extending the table and indicating in ^rtiich COMMON the new 
variable resides. Also include the address within that COMMON. A study 
of Figure 12 will indicate the simple coding change required to add a new 
labeled COMMON to the table. 
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8.3 Adding QuAOtities to Output Set 

The general mechanism by which the output signal set is defined has been 
discussed elsewhere in this volume as well as in Volume III. The key item 
to remember is that anj' program variable defined in blank CCM^!0N can be chosen 
for the signal set. 

Blank COMMON contains all the program variables considered output 
material as of this writing, A mechanism for extending the set has been 
provided, namely, the vector OUTD at the end of CCMMOS and the subprogram 
AUXOUT. 

Ihe auxiliary output routine, AUXOUT, defines eleaents of the output 
vector, OUTD. To include a new output, define a heretofore undefined element 
of OUTD. Ihus, its output identification is defined. It should be stressed 
that blank COMMON should not be redefined in any way. Any alteration of 
the order of variables in COM<ON will invalidate current output identifications. 



9. PROGRAM CROSS REFERENCE MAP 


The collection of subprograms »rtiich constitute REXOR II ccaatain many 
cross-referenced COMMON blocks and subroutines. Figujre 30 is presented as an 
aid to developing the source smd usage of any particular item. The vertical 
listing in this figure gives all the routines, and names or unnamed (one) 
CO'ilON blocks used in REXOR II. A second column, headed 'T* , gives the use 
of the entry. The coding used is *M' for main program, *S* for subroutine, 

'E' for entry of a subroutine, 'C* for common block. The list of names is 
in alphabetical order and numbered. The alphabetizing is by main to sub- 
routines with a subalphabetizing of entry points under each subroutine. The 
COMMON blocks are listed last except for subroutines not included in the 
source deck, "niese subroutines are usually part of the computer operating 
(precompiled) package, and not particularly associated with REXOR II. Ilie 
numbering is repeated horizontally, and corresponds to the vertical name list. 

Ihe vertical list on the left-hand side is the csilling or active routine 
or element, and the horizonted line lists the routines called or referenced. 
Numbers at grid intersections show there is a reference and the level ~f 
reference. One indicates a direct reference. Two or three show there are 
one or two intermediate references, respectively. Bote that a subroutine name 
will show all the references to all the entry points boxinded by that name. 

Sy elimination, the references associated with the subroutine name up to the 
point of the first entry name can be determined. 
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Figure 30. - Connnon eubroutin* directory. 
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Figure 30. - Concluded 



10. PROGRAM LISTIHG 


Due to the large number of pages in the REXOR II program source listing, 
this material is handled under a separate binder. Copies may be obtained from 
the distributing agency, RASA, Langley Research Center, Hampton, Virginia. 
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